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and  Their  Applications 
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Madison,  Wisconsin 
June  25  -  July  1,  1989 


The  Eighteenth  Conference  on  Stochastic  Processes  and  Their  Applications  was 
held  at  the  University  of  Wisconsin  -  Madison,  June  25  -  July  1,  1989.  The 
Conference  included  140  participants  from  20  countries.  One-hour  invited 
lectures  were  given  by  the  following: 


K.  Alexander  (USA) 

V.  Anantharam  (USA) 

P.  Baxendale  (USA) 

K.  Burdzy  (USA) 

M.  Chaleyat-Maurel  (France) 
R.  Dobrushin  (USSR) 

J.  Le  Gall  (France) 

G.  Grimmett  (UK) 


B.  Hajek  (USA) 

P.  Hsu  (USA) 

P.  March  (Canada) 

G.  Papanicolaou  (USA) 
R.  Schonmann  (Brazil) 
M.  Talagrand  (France) 
K.  Uchiyama  (Japan) 


In  addition,  approximately  70  contributed  lectures  were  presented.  Abstracts 
of  the  invited  and  contributed  lectures  will  be  published  in  Stochastic 
Processes  and  Their  Applications. 

The  program  conmittee  consisted  of  D.  Griffeath  (chairman) ,  A.  de  Acosta,  M. 
Bramson,  P. Glynn,  F.  Kelly,  M.  Pinsky,  and  R.  Williams.  Local  arrangements  were 
handled  by  D.  Griffeath,  T.  Kurtz,  J.  Kuelbs  and  V.  Leatherberry.  The 
conference  was  organized  under  the  auspices  of  the  Committee  for  Conferences  on 
Stochastic  Processes  and  the  Bernoulli  Society,  with  sponsorship  and  support 
from  the  Department  of  Mathematics,  University  of  Wisconsin  -  Madison. 

Financial  support  was  provided  by  the  Air  Force  Office  of  Scientific  Research, 
the  Army  Research  Office,  the  National  Security  Agency,  the  Office  of  Naval 
Research,  and  the  National  Science  Foundation.  This  financial  support  provided 
travel  and  subsistence  for  the  invited  speakers  and  20  graduate  students  and 
young  faculty  who  did  not  have  other  research  support. 
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B.  Hajek  Load  balancing  for  a  random  graph 

Univ.  of  Illinois 
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Abstracts  of  Invited  Lectures 


LARGE  FINITE  CLUSTERS  IN  TWO-DIMENSIONAL  PERCOLATION 

Kenneth  S.  Alexander 

Department  of  Mathematics 
University  of  Southern  California 

Jennifer  T.  Chayes 
Department  of  Mathematics 
UCLA 

Lincoln  Chayes 
Department  of  Mathematics 
UCLA 

For  independent  bond  percolation  in  two  dimensions  at  density  p  above  the  critical 

point,  some  questions  relating  to  large  Finite  clusters  are  considered.  The  probabilities 

Pj  [0  «  x]  of  connection  by  a  path  of  occupied  bonds  at  density  1  -  p  are  used  in 

2 

constructing  a  natural  norm  gp(  • )  on  R  .  It  is  shown  that  the  probability  Pp[  |  C(0)  |  =  N] 
that  the  cluster  of  the  origin  has  finite  size  N  behaves  like  exp(-R(p)N^)  as  N-», 
where  R(p)  involves  the  correlation  length,  the  percolation  probability,  and  an  isoperimetric 
constant:  the  minimum  gp-length  of  any  curve  enclosing  unit  area.  If  W(p)  is  the  region 
enclosed  by  the  corresponding  shortest  curve,  then  conditionally  on  [  |C(0)  |  =  N],  the  cluster 
C(0)  of  the  origin  closely  approximates  W(p)  in  shape,  with  probability  approaching  one  as 
N  - 
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THE  TRANSIENT  BEHAVIOR  OF  NETWORKS  OF  QUEUES 


V.  Anantharam 

School  of  Electrical  Engineering 
313  Phillips  Hall 
Cornell  University 
Ithaca,  New  York  14853 

A  closed  Jackson  network  is  a  network  of  •  /M/1  queues,  where  service  at  each  of  the 
nodes  is  FCFS  and  customers  leaving  a  node  are  routed  to  one  of  the  other  nodes  according 
to  a  stochastic  matrix,  the  routing  decisions  being  independent  of  the  service  times.  We  will 
discuss  the  probabilistic  evolution  of  such  a  network  when  there  is  a  large  number  of 
customers  circulating.  In  particular,  we  will  prove  the  existence  of  a  threshold  phenomenon  in 
the  approach  to  stationarity.  The  estimates  going  into  the  proof  also  allow  one  to  derive  a 
useful  rule  of  thumb  for 
problem  of  considerable 


the  buffer  allocation  problem  for  open  networks  of  queues  —  a 
applied  interest.  This  derivation  will  also  be  presented. 
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STATISTICAL  EQUILIBRIUM  AND  MIXING  FOR  STOCHASTIC  FLOWS 

OF  DIFFEOMORPHISMS 

Peter  H.  Baxendale 
University  of  Southern  California 

Let  {^t :  £0}  be  a  stochastic  flow  of  diffeomorphisms  of  a  compact  manifold  M. 
That  is,  :  £0}  is  a  Diff(M)  valued  process  with  continuous  sample  paths,  started  at  Id, 
with  stationary  left  increments  (for  (Ks^t<«).  Let  P(M)  denote  the  space  of  Borel 

probability  measures  on  M.  Then  pt  =  p£^  (peP(M))  induces  a  Markov  process  {pt :  0} 

in  P(M). 

If  p  is  the  unit  mass  at  x  in  M  then  {pt :  t£0)  is  naturally  jomorphic  to  the 
one-point  motion  {^(x) :  £0}  in  M.  More  interestingly,  if  p  is  stationary  for  the 
one-point  motion  then  (using  martingale  convergence)  p{  converges  in  distribution  as  t-*« 
to  the  so-called  statistical  equilibrium.  Furthermore,  in  a  large  class  of  examples  the 
statistical  equilibrium  consists  almost  surely  of  measures  of  Hausdorff  dimension  strictly 
between  0  and  dim  (M).  [Results  due  to  LeJan  and  Ledrappier  and  Young.] 

In  this  talk  we  study  the  behavior  as  t  -*  «  of  p£  for  arbitrary  p  in  P(M).  In 
particular  we  investigate  when  (and  at  what  rate)  Pj  t  -  ^  t  converges  to  zero  as  t  -* »  for 
distinct  p^  and  p^  in  P(M). 
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NON-INTERSECTION  EXPONENTS  FOR  BROWNIAN  PATHS 

Krzysztof  Burdzy 
University  of  Washington 

Let  X  and  Y  be  independent  n-dimensional  Brownian  motions,  n  £  2,  X(0)  =  0, 

|  Y (0)  |  =  e,  and  let  pn(E)  =  P(X[0,1]  n  Y[0,1]  =  0).  Asymptotic  estimates  (when  e  -♦  0)  of 
pn(£)  and  some  related  probabilities  will  be  discussed.  Applications  to  geometry  of  Brownian 
fractals  will  be  given. 
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SAMPLE  PATH  PROPERTIES  OF  A  CLASS  OF  HYPOELLIPTIC 
DIFFUSION  PROCESSES 

Mireille  Chaleyat— Maurel 
University  Paris  VI 

3 

We  study  a  simple  class  of  hypoelliptic  non-elliptic  diffusion  processes  in  R  .  We 
give  precise  estimates  for  the  Green  function  of  the  process  and  the  capacity  of  small  compact 
sets.  The  proofs  rely  on  a  study  of  the  local  behavior  of  the  process,  developed  along  the  lines 
ofAzen  and  on  some  bounds  for  the  Greei  iction  which  have  bee  ’,  recently  derived  by 
analytic  methods.  Our  estimates  can  be  applies  various  sample  path  iperties,  such  as  the 
limiting  behavior  of  the  volume  of  a  small  tubuior  neighborhood  of  the  path,  or  the  existence 
of  multiple  points.  We  emphasize  the  differences  with  the  elliptic  case.  In  particular,  the 
volume  of  a  tubular  neighborhood  of  the  path  is  shown  to  be  smaller  than  for  an  elliptic 
diffusion. 
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GIBBS  SHAPE  OF  A  DROPLET 


Roland  Dobrushin 

Institute  of  Problems  in  Information  Transmission 

We  consider  the  Gibbs  distribution  for  a  low-temperature  two-dimensional  Ising 
model  with  a  periodic  boundary  condition  in  volume  V.  We  introduce  an  additional  condition 
that  the  total  spin  takes  a  fixed  value  M(V).  Supposing  that  limM(V)/|V|  =m  if  |V| 
goes  to  infinity  where  |  m  j  >  1/2,  we  consider  the  asymptotic  limit  of  the  corresponding 
conditional  distribution.  It  turns  out  that  it  converges  to  the  situation  when  there  is  a  droplet 
of  one  phase  inside  the  other  phase.  This  limit  droplet  has  a  fixed  shape  which  is  called  the 
Wulff  shape  in  statistical  mechanics.  This  result  was  obtained  by  R.  Kotezky,  S.  Shlosman 
and  myself.  In  the  talk  many  open  questions  connected  with  the  generalization  of  this  result 
will  also  be  discussed. 
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SLAB  RESULTS  FOR  PERCOLATION  AND  CONTACT  PROCESSES 


Geoffrey  Grimmett 
University  of  Bristol 

Of  the  many  major  open  problems  for  infinite  particle  systems  and  random  media, 
some  the  the  most  accessible  concern  the  behavior  of  processes  at  their  critical  points.  The 
critical  branching  process  dies  out  (barring  the  special  case),  but  is  the  same  true  for 
percolation  and/or  the  contact  process,  say?  If  this  is  so,  then  what  about  critical  exponents...? 

Related  he  question  about  behavior  at  the  critical  point  is  the  question  of  deciding 
to  what  degree  a  process  in  d  dimensions  may  be  'approximated'  by  a  process  living  on  an 
essentially  two-dimensional  subset  (there  are  special  techniques  for  processes  in  two 
dimensions). 

We  arrive  thus  at  the  following  matters: 

PERCOLATION:  Is  there  percolation  when  the  density  of  open  bonds  equals  its  critical  value 

p  ?  Is  the  limit  of  the  critical  probabilities  of  slabs  equal  to  p  ? 
c  c 

CONTACT  PROCESS:  Does  the  contact  process  survive  when  the  infection  rate  equals  its 

deal  value  X  ?  What  does  the  infected  region  in  space-time  look  like,  for  general  starting 
c 

sets  and  infection  rates? 

Recent  process  in  slab  technology  has  led  to  sensible  answers  to  almost  all  of  these 
questions,  and  hence  resolutions  of  several  problems  of  interest  for  supercritical  percolation 
and  contact  processes. 
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LOAD  BALANCING  FOR  A  RANDOM  GRAPH 


B.  Hajek* 

University  of  Illinois  —  Urbana 

A  set  of  M  resource  locations  and  a  set  of  cxM  consumers  are  given.  Each  consumer 
requires  a  specified  amount  of  resource  and  is  constrained  to  obtain  the  resource  from  a 
specified  subset  of  locations.  The  problem  of  assigning  consumers  to  resources  so  as  to 
balance  the  load  among  the  resource  locations  as  much  as  possible  is  considered.  It  is  shown 
that  there  are  assignments,  termed  uniformly  most  balanced  assignments,  that  simultaneously 
minimize  certain  symmetric,  separable,  convex  cost  functions. 

The  distribution  function  of  the  load  at  a  given  location  for  a  uniformly  most  balanced 
assignment  is  studied  assuming  that  the  set  of  locations  each  consumer  can  use  is  random.  An 
asymptotic  lower  bound  (conjectured  to  be  asymptotically  exact)  on  the  distribution  function  is 
given  for  M  tending  to  infinity,  and  an  upper  bound  is  given  on  the  probable  maximum  load. 

It  is  shown  that  there  is  typically  a  large  set  of  resource  locations  that  all  have  the 
maximum  load,  and  that  for  large  average  loads  the  maximum  load  is  near  the  average  load. 

* 

Sponsored  by  the  Joint  Services  Electronics  Program,  Contract  No.  N00014-84-C-0149. 
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PROBABILISTIC  METHODS  IN  DIFFERENTIAL  GEOMETRY 


Pei  Hsu 

Northwestern  University 

We  will  introduce  briefly  Riemannian  Brownian  motion  and  the  minimal  heat  kemal  on 
a  Riemannian  manifold  and  discuss  how  the  Ricci  curvature  and  sectional  curvature  affects  the 
behavior  of  Brownian  motion.  We  then  show  how  Brownian  motion  can  be  used  to  prove 
some  old  and  new  results  in  differential  geometry,  including  the  existence  of  harmonic 
functions  on  Cartan— Hadamard  manifolds,  small  time  asymptotic  behavior  of  the  heat  kemal, 
the  stochastic  completeness  and  Feller  property  of  the  heat  semigroup,  and  heat  kernel  on 
noncomplete  Riemannian  manifolds. 


14 


A  TRAJECTORIAL  CONSTRUCTION  OF  SOME  MEASURE-VALUED  PROCESSES 

Jean— Francois  Lc  Gall 
University  Paris  VI 

We  give  a  trajectorial  construction  of  some  measure-valued  Markov  processes  which 
have  been  studied  recently  by  Dawson,  Dynkin,  Perkins  and  others.  These  processes,  called 
superprocesses  by  Dynkin,  were  originally  constructed  as  weak  limits  of  systems  of  branching 
diffusions.  Our  approach  uses  the  theory  of  excursions  of  linear  Brownian  motion,  and  the 
notion  of  the  branching  tree  associated  with  a  Brownian  excursion.  It  yields  a  simple 
probabilistic  description  of  the  nonlinear  semi-group  that  appears  in  the  study  of  the 
measure-valued  process.  Another  advantage  of  our  method  is  as  follows.  If  we  interpret  the 
value  of  the  process  at  a  given  time  t  as  a  measure  on  a  cloud  of  particles,  our  approach 
makes  it  possible  to  consider  the  induced  measure  on  the  set  of  paths,  stopped  at  time  t,  of 
these  particles.  Some  questions  related  to  the  genealogy  of  the  particles  can  be  treated  in  this 
way,  and  one  may  expect  applications  to  sample  path  properties  of  the  type  investigated  by 
Perkins. 
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REMARKS  ON  THE  WTTTEN-SANDER  MODEL 


Peter  March 
Ohio  State  University 

The  Witten— Sander  model,  or  diffusion  limited  aggregation,  is  a  simply  stated  cluster 
growth  model  on  the  plane  lattice.  We  make  an  indirect  statement  about  the  growth  of  a 
scaled  sequence  of  clusters  of  a  similar  type  by  a  comparison  with  the  ill-posed  Hele-Shaw 
problem  of  fluid  mechanics. 
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DIRECT  AND  INVERSE  PROBLEMS  FOR  PULSE  REFLECTION 
BY  RANDOMLY  LAYERED  MEDIA 

George  Papanicolaou 

Courant  Institute 
251  Mercer  Street 
New  York,  NY  10012 


Consider  a  pulsed  acoustic  or  elastic  wave  impinging  on  a  material  that  is 
inhomogeneous,  the  earth's  crust  for  example.  The  reflected  signal  is  measured  at  the  surface 
and  it  looks  very  noisy,  like  a  seismogram.  What  information  about  the  medium  can  we 
extract  from  this  signal?  In  reflection  seismology  this  is  a  classic  question  that  has  received  a 
great  deal  of  attention  over  the  last  forty  years.  It  has  also  received  a  lot  of  attention  as  a 
difficult  problem  in  mathematics:  the  inverse  scattering  problem.  We  have  begun  a  study  of 
this  question  as  a  problem  in  stochastic  processes.  We  model  the  reflecting  material  as  a 
random  medium,  characterize  all  possible  reflected  signals  as  a  particular  class  of 
nonstationary  stochastic  processes  and  then  pose  the  recovery  of  material  properties  as  a 
statistical  estimation  problem  for  this  class  of  processes.  We  have  also  conducted  extensive 
numerical  simulations  to  test  the  theory;  the  estimation  algorithms  in  particular. 
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ON  THE  BEHAVIOR  OF  SOME  CELLULAR  AUTOMATA 

Roberto  H.  Schonmann 

Sao  Paulo  University 
Brazil 

We  study  the  behavior  of  some  cellular  automata  related  to  Bootstrap  Percolation. 

Each  site  of  the  Lattice  2^  is  in  state  0  or  1.  Sites  which  are  in  state  1  remain  so  forever, 
while  sites  that  are  in  state  0  may  change  to  state  1  deterministically  according  to  the  states 
of  their  neighbors,  following  some  specified  rules.  All  the  randomness  is  in  the  initial  state. 

A  basic  exam,  is  the  one  :  lich  at  each  integer  time  a  site  which  is  in  state  0 
flips  to  state  1  if  at  least  half  of  its  jighbors  are  in  state  1 .  We  will  show  that  in  any 
dimension  and  for  any  initial  i.i.d.  distribution  with  positive  density  every  site  will  eventually 
turn  to  state  1.  We  will  then  discuss  how  fast  this  happens.  The  methods  used  will  provide 
some  simple  examples  of  techniques  such  as  Peierls  estimates  and  renormalization. 
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APPLICATIONS  OF  ISOPERIMETRIC  INEQUALITIES  TO  SUMS 
OF  INDEPENDENT  RANDOM  VARIABLE. 

M.  Talagrand 
University  Paris  VI 

We  introduce  a  new  method  to  estimate  the  tails  of  sums  of  independent  Banach  space 
valued  random  variables.  The  method  relies  on  two  isoperimetric  inequalities  of  independent 
interest,  and  yields  sharp  results. 


ON  DERIVATION  OF  THE  BOLTZMANN  EQUATION 
FROM  A  DETERMINISTIC  MOTION  OF  MANY  PARTICLES 

K.  Uchiyama 

Hiroshima  University 
Hiroshima,  Japan 

We  introduce  a  dynamical  system  of  "hard"  squares  (having  velocities  from  a  finite  set) 
that  seemingly  corresponds  to  the  Broadwell  model  of  the  B.  eq.  (Boltzmann  equation)  and 
discuss  the  Grad  limit  for  it.  For  the  hard-sphere  model  O.  E.  Lanford  derived  the  B.  eq.  from 
the  BBGKY  hierarchy  for  short  times  by  taking  the  Grad  limit.  For  the  present  model  an 
analogy  is  observed  at  the  formal  level:  the  first  equation  of  the  BBGKY  hierarchy  is  reduced 
to  the  corresponding  B.  eq.  (Broadwell  model)  by  formally  passing  to  the  Grad  limit  The 
actual  limit  which  can  be  shown  to  exist  for  short  times,  however  does  not  solve  the  B.  eq. 

It  is  commonly  understood  that  the  B.  eq.  is  valid  under  the  hypothesis  of  Molecular 
Chaos,  but  both  the  precise  meaning  of  the  hypothesis  and  the  way  how  it  assures  tye  B.  eq. 
remain  to  be  clarified.  H.  Grad  pointed  out  that  the  chaos  property  for  two  particles  is 
destroyed  for  post-collisional  configurations  and  in  order  to  assure  the  validity  of  the  B.  eq.  it 
is  needed  only  for  those  configurations  in  which  particles  are  about  to  collide.  Now  is  it 
plausible  to  suppose  that  if  recollisions  are  negligible  so  that  two  particles  that  are  about  to 
collide  with  each  other  have  rarely  made  interaction  through  their  "ancestorial"  particles,  then 
they  must  be  uncorrelated?  Our  model  serves  as  a  counter  example  to  this  supposition:  the 
recollisions  are  negligible,  still  such  two  particles  are  correlated. 
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Contributed  Papers 
Monday,  June  26 


Track  A 
RoomB207 

Queueing  Models  and  Related  Topics 


2:00  -  2:20  p.m. 

R.  J.  Williams 

On  the  quasireversibility  of  a 
multidass  Brownian  service  station 

2:30  -  2:50  p.m. 

M.  L  Reiman 

Light  and  heavy  traffic  limits  for  two 
finite  capacity  queueing  systems 

3:00  -  3:20  p.m. 

E.  V.  Krichagina 

The  heavy  traffic  approximation  for 
closed  exponential  queueing  networks 

4:00  -  4:20  p.m. 

J.  Dshalalow 

A  duality  principle  in  doubly  controlled 
processes  of  servicing  machines 
with  reserve  replacement 

4:30  -  4:50  p.m. 

D.  Perry 

On  arrivals  and  departures  that  see 
the  same  average  in  bulk  queues 

Track  B 
Room  B215 
Limit  Theorems 


2:00  -  2:20  p.m. 

M.  S.  Bingham 

An  approximate  martingale  central 
limit  theorem  in  Hilbert  space 

2:30  -  2:50  p.m. 

M.  Peligrad 

On  the  tightness  of  the  partial  sums 
of  a  O- mixing  sequence 

3:00  -  3:20  p.m. 

W.  Vervaat 

Loglog  laws  for  geometric  and 
full  sequences 

4:00  -  4:20  p.m. 

W.  Bryc 

Limit  theorems  via  large  deviations 

4:30  -  4:50  p.m. 

V.  H.  de  la  Pena 

Precision  bounds  on  the  Lp  norm  of  a 
quadratic  form  of  independent 
mean  zero  random  variables 
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Track  C 
Room  B223 

Biological  and  Physical  Science  Modeling 


2:00  -  2:20  p.m. 

S.  N.  Ethier 

The  distribution  of  the  frequencies  of 
age-ordered  alleles  in  a  diffusion  model 

2:30  -  2:50  p.m. 

P.  Taylor 

Product  form  in  migration  processes 
with  batch  movements 

3:00  -  3:20  p.m. 

D.  Blount 

A  law  of  large  numbers  and  central 
limit  theorem  for  a  stochastic  model  of  a 
nonlinear  chemical  reaction  with  diffusion 

4:00  -  4:20  p.m. 

F.  Ball 

Aggregated  Markov  processes  with 
negative  exponential  time  interval  omission 

4:30  -  4:50  p.m. 

G.  J.  Morrow 

Large  deviation  results  for  a  class  of 

Banach  space  valued  Markov  chains  with 
application  to  population  genetics 

Track  D 

Room  B239 

Diffusion  and  Stochastic  Differential  Equations 

2:00  -  2:20  p.m. 

J.  Picard 

Non-linear  filtering  with  small 
noise  coefficients 

2:30  -  2:50  p.m. 

J.-L.  Jiang 

On  the  exchangeability  of  differentiation 
and  integration  with  applications  to 
stochastic  p.d.e.'s 

3:00  -  3:20  p.m. 

M.  Pinsky 

Central  limit  and  recurrent  behavior 
of  diffusions  in  a  random  environment 

4:00  -  4:20  p.m. 

X.  Mao 

Exponential  stability  for  stochastic 
differential  equations  with  respect  to 
semimartingales 

4:30  -  4:50  p.m. 

G.  S.  Katzenberger 

Semimartingales  forced  onto  a  manifold 
by  a  large  drift 

Contributed  Papers 
Tuesday,  June  27 


Track  A 
Room  B207 

Topics  in  Stochastic  Processes 


2:00  -  2:20  p.m. 
2:30  -  2:50  p.m. 

3:00  -  3:20  p.m. 

4:00  -  4:20  p.m. 
4:30  -  4:50  p.m. 


S.  Cambanis  Trapezoidal  Monte  Carlo  integration 

H.  L.  Hurd  Correlation  theory  of  almost 

periodically  correlated  processes 

L.  D.  Pitt  A  probabilistic  analysis  of  Lipschitz  points 

for  non-differentiable  A*  functions 

(To  be  announced) 

(To  be  announced) 


Track  B 
Room  B215 

Topics  in  Extremal  Theory 


2:00  -  2:20  p.m. 
2:30  -  2:50  p.m. 

3:00  -  3:20  p.m. 
4:00  -  4:20  p.m. 


G.L.  O'Brien 

J.  M.  P.  Alb  in 

K.  F.  Turkman 
M.  G.  Hahn 


Loud  shot  noise 

On  extremal  theory  for  differentiable 
stationary  processes 

Sums  and  maxima  of  stationary  sequences 

Max  infinitely  divisible  and  max  stable 

sample  continuous  processes 


4:30  -  4:50  p.m.  D.  Ross 


A  representation-independent  integral 
for  set-valued  random  variables 
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Track  C 
Room  B223 

Estimation  and  Statistical  Inference 


2:00  -  2:20  p.m. 

A.  G.  Gadziev 

Regression  models  with  increasing 
number  of  unknown  parameters 

2:30-  2:50  p.m. 

K.  J.  Utikal 

Non  parametric  inference  for  a 
doubly  stochastic  Poisson  process 

3:00  -  3:20  p.m. 

A.  K.  Basu 

Local  asymptotic  mixed  normality  of 
log-likelihood  based  on  stopping  times 
and  its  rate  of  convergence 

4:00  -  4:20  p.m. 

M.  A.  A.  Turkman 

Prediction  for  an  autoregressive  process 
with  exponential  white  noise 

4:30  -  4:50  p.m. 

N.  Ebrahimi 

T.  Ram  al  ingam 

Estimation  of  reliability  of  stress-strength 
models  based  on  sample  paths 

Track  D 
Room  B239 

Diffusion  and  Partial  Differential  Equations 


2:00  -  2:20  p.m. 

2:30  -  2:50  p.m. 
3:00  -  3:20  p.m. 
4:00  -  4:20  p.m. 


G.  Galeazzi 

C  Graham 

G.  Giroux 

A.  Korzeniowski 


Uniqueness  of  viscosity  solutions  of 
Hamil .  .  ^Jacobi-Bellman  equations  for 
controlled  diffusions  on  finite-dimensional 
Riemannian  manifolds  with  boundary 

Nonlinear  limit  for  a  system  of  diffusing 
particles  which  alternate  between  two  states 

A  geometric  rate  of  convergence  for  the 
Maxwell  gas  with  angular  cut-off 

A  viscosity  algorithm  for  the  scalar 
conservation  law 


4:30  -  4:50  p.m.  S.  Roy 


Stochastic  geometry  and  relativistic 
path  integrals 
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Contributed  Papers 
Thursday,  June  29 


Track  A 
Room  B207 

Branching  Processes  and  Point  Processes 


2:00  -  2:20  p.m. 

P.  Broberg 

2:30  -  2:50  p.m. 

A  Chen 

3:00  -  3:20  p.m. 

P.  Blackwell 

4:00  -  4:20  p.m. 

G.  Nieuwenhuis 

4:30  -  4:50  p.m. 

A.  Bose 

Dependencies  in  branching  processes 

Markov  branching  processes  with 
instantaneous  immigration 

Stochastic  models  for  the  evolution  of 
social  behaviour  based  on  the  resource 
dispersion  hypothesis 

Asymptotics  for  point  processes 

A  fluctuations  limit  for  scaled  age 
distributions  and  weighted  Sobolev  spaces 


Track  B 
Room  B215 

Set-indexed  Processes  and  Levy  Processes 


2:00  -  2:20  p.m. 

M.  Ossiander 

The  size  of  domains  of  set-indexed 

L6vy  processes 

2:30  -  2:50  p.m. 

J.  Mijnheer 

U-statistics  and  the  double  stable  integral 

3:00  -  3:20  p.m. 

J.  E.  Yukich 

Some  new  Vapnik-Chervonenkis 
classes  of  sets 

4:00  -  4:20  p.m. 

J.  Levy 

Measuring  the  long-term  dependence  of 
some  stationary  stable  processes 

4:30  -  4:50  p.m. 

R.  A.  Doney 

Hitting  probabilities  for  spectrally 
positive  Levy  processes 
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Track  C 
Room  B223 

Random  Graphs  and  Probabilistic  Algorithms 


2:00  -  2:20  p.m. 

J.  van  den  Berg 

Some  probabilistic  problems  arising 
from  an  optimal-path  algorithm 

2:30  -  2:50  p.m. 

M.  Cocan 

Upon  a  possibility  of  generation  of 
random  graphs 

3:00  -  3:20  p.m. 

F.  Avram 

Convergence  to  the  normal  distribution 
by  graph  methods 

4:00  -  4:20  p.m. 

(To  be  announced) 

4:30  -  4:50  p.m. 

(To  be  announced) 

Track  D 

Room  B239 

Brownian  Motion  and  Related  Topics 

2:00  -  2:20  p.m. 

L.  C  G.  Rogers 

Brownian  motion  in  a  wedge  with 
variable  skew  reflection 

2:30  -  2:50  p.m. 

J.  van  der  Weide 

Random  walk  approximations  of 
skew  and  sticky  Brownian  motion 

3:00  -  3:20  p.m. 

R.  J.  Vanderbei 

New  results  for  optimal  switching 
between  a  pair  of  Brownian  motions 

4:00  -  4:20  p.m. 

M.  A,  Amir 

Brownian  motion  as  the  strong  limit  of  the 
difference  of  two  Poisson  processes  with 
an  application  to  M/M/1  queues 

4:30  -  4:50  p.m. 

T.  Lindstrom 

Brownian  motion  on  nested  fractals 
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Contributed  Papers 
Friday,  June  30 


Track  A 
Room  B207 

Random  Times  and  Related  Topics 


2:00  -  2:20  p.m. 
2:30  -  2:50  p.m. 

3:00  -  3:20  p.m. 
4:00  -  4:20  p.m. 

4:30  -  4:50  p.m. 


S.  Krone 
O.  Kallenberg 

C  Ching-Sung 
G.  Roberts 

(To  be  announced) 


Local  times  for  superdiffusions 

Random  time  change  and  int  il 
representations  for  marked  s;  ping  times 

On  some  inequalities  of  local  time 

Approximations  for  diffusion  process 
boundary  hitting  times 


Track  B 
Room  B215 

Topics  in  Probabilistic  Modeling 


2:00  -  2:20  p.m. 
2:30  -  2:50  p.m. 

3:00  -  3:20  p.m. 


W.  A.  Massey 
W.  Willinger 

P.  R.  Vittal 


Steady  state  analysis  for  a  series  of 
queueing  networks  with  blocking 

Equivalent  martingale  measures  and 
no-arbitrage  in  stochastic  securities 
market  models 

Expected  amount  of  overflow  in 
a  finite  dam 


4:00  -  4:20  p.m.  V.  T.  Stef anov 

4:30  -  4:50  p.m.  G.  V.  Orman 


On  the  finiteness  of  the  moments  of 
first-crossing  times  from  nonlinear 
renewal  theory 

A  problem  of  expectation  and  variance 
for  a  class  of  formal  languages 
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Track  C 
Room  6223 

Particle  Systems  and  Percolation 


2:00  -  2:20  p.m. 

R.  Fisch 

Cyclic  cellular  automata 
and  related  processes 

2:30  -  2:50  p.m. 

D.  Chen 

Finite  reversible  nearest-partide 
systems  in  inhomogeneous  and 
random  environments 

3:00  -  3:20  p.m. 

N.  Madras 

Using  renewal  theory  to  investigate 
self-avoiding  walks 

4:00  -  4:20  p.m. 

J.  C  Wierman 

The  substitution  method  for  bond 
percolation  critical  probability  bounds 

4:30  -  4:50  p.m. 

(To  be  announced) 
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Abstracts  of  Contributed  Papers 


ON  EXTREMAL  THEORY  FOR  DIFFERENTIABLE  STATIONARY  PROCESSES 

J.  M.  P.  Albin 

University  of  Lund 
Sweden 

Let  {X(t)}^Q  be  a  stochastically  differentiable  strictly  stationary  stochastic  process, 
and  let  be  a  decreasing  family  of  rare  subsets  of  IRm,  i.e.  P{X(0)€Au)  -»  0,  as 

u  -*  «.  We  give  a  general  method  to  derive  an  asymptotic  expression  for  the  probability 
P{X(t)  e  Au>  for  some  (Kt^h)  as  u  -*  «.  In  particular,  letting  Au  =  (u,«0  and  m  =  1,  the 
method  suggests  a  way  to  study  the  "usual"  extremal  problem. 

Our  approach,  leaning  more  directly  on  the  differentials  of  X(t)  than  through  their 
impact  on  the  moments  of  the  number  of  upcrossings,  yields  a  clustering  allowing  extremal 
theory  with  conditions  formulated  in  terms  of  f.d.d.'s  of  order  three  and  lower,  while  the 
traditional  upcrossing  approach,  which  is  not  clustering  allowing,  needs  conditions  involving 
f.d.d.'s  of  order  four. 
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PREDICTION  FOR  AN  AUTOREGRESSIVE  PROCESS 
WITH  EXPONENTIAL  WHITE  NOISE 


Maria  AnttSnia  Amaral  Turkman 

CEA  and  Dep.  Estatfstica 
University  of  Lisbon 

Autoregressive  processes  with  non-Gaussian  white  noise  are  useful  for  modelling  a 
wide  range  of  phenomena  which  do  not  allow  negative  values.  Such  processes  were  studied 
by  Gaver  and  Lewis  (1980),  Bell  and  Smith  (1986)  among  others.  Andel  (1988)  considers  the 
autoregressive  model  Xt  =  pYt_^  +  Y£  where  0  £  p  <  1,  Y{  are  independent  exponential 
random  variables  with  mean  value  9~*  and  Xj  is  exponential  with  mean  value  (0(1— p))~*. 
This  model  was  studied  from  a  Bayesian  point  of  view  by  Amaral  Turkman  (1989).  In  this 
paper  we  address  the  problem  of  predicting  certain  type  of  events  which  may  be  of  interest  in 
a  practical  situation  modeled  by  this  process. 

References 

Amaral  Turkman,  M.  A.  (1989)  Bayesian  Analysis  of  an  autoregressive  process  with 
exponential  white  noise.  Submitted  to  Statistics. 

Andel,  J.  (1988)  On  AR(1)  processes  with  exponential  white  noise.  Commun.  Statist.  Theory 
Meth.  17  (5),  1481-1495. 

Bell,  C.  B.  &  Smith,  E.  P.  (1986)  Inference  for  non-negative  autoregressive  schemes. 
Commun.  Statist.  Theory  Meth.  15  (8),  2267-2293. 

Gaver,  D.  P.  &  Lewis,  P.  A.  (1980)  First  order  autoregressive  gamma  sequences  and  point 
processes.  Adv.  Appl.  Prob.  12,  727-745. 
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BROWNIAN  MOTION  AS  THE  STRONG  LIMIT  OF  THE  DIFFERENCE  OF  TWO 
POISSON  PROCESSES  WITH  AN  APPLICATION  TO  M/M/1  QUEUES 

Madjid  A.  Amir  and  Frank  B.  Knight 
University  of  Illinois 


We  consider  two  independent  Poisson  processes  An(t)  and  Bn(t)  with  rates  Xn  and 

An(t)-B(t) 

|i  respectively  defined  on  (QJvJPJ,  for  n  =  1,2,3,... .  Let  R_(t)  = - - — .  R  (t)  is 

n  n  n  n  n  2**  ^ 


then  a  birth  and  death  process  with  state  space  the  lattice  {— Jc  =  ±l^h2r±3,...},  and  birth  rate 

2n 

X  u 

(death  rate  resp.)  i  ;rv  (vn  +LL  ~»  resp.).  We  show  the  existence  of  a  probability  space 
'Si  ^n  A,n+^n 

(O^.F^.P^)  (projective  limit  space)  such  that  =  {co  =  ((DpC^*0^*-*^.- )»®n  e  &n)  and 

Poo(Sn)  =  Pn(Sn)  for  Sn  e  Fn  and  n  =  1,2,3 .  We  show  that 


P  jlim  sup 


An(t)-Bn(t) 


n-*«  te[0, 1]  2n 


-  W(t)  |  =  0|  =  1 


where  W(t)  is  standard  Brownian  motion  (Brownian  with  drift  -ln(a),  resp.)  when 

=  p-n  =  2z  (Xn=2  and  =  2  "a1  ,  resp.),  where  a  is  any  positive  real  number. 

We  thus  construct  a  Brownian  motion  with  and  without  drift  as  a  strong  limit  of  a  sequence  of 
birth  and  death  processes.  Note  that  the  weak  convergence  in  the  space  D  will  follow. 

As  many  processes  arising  in  applied  probability  can  be  modeled  as  functionals  of  the 
difference  of  two  Poisson  processes,  this  result  has  applications  in  areas  such  as  queueing 
theory  and  storage  theory.  A  natural  application  is  to  M/M/1  queues:  Let  Qn(t)  be  the 
queue  length  process  of  an  M/M/1  queue  with  either  set  of  rates  above,  then: 


1 i m  sup  J 
n-*»  te[0, 1] 


Q„(0 


-  W(t)  I  =  0=1 


where  W(t)  =  W(t)  -  infs^tW(s)  is  reflected  Brownian  motion.  This  is  a  much  stronger  result 
than  the  usual  weak  convergence.  Some  rates  of  convergence  have  been  derived. 
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CONVERGENCE  TO  THE  NORMAL  DISTRIBUTION  BY  GRAPH  METHODS 


Florin  Avram 
Northeastern  University 

When  establishing  convergence  to  the  normal  distribution  of  sums  of  Wick  products  of 
a  stationary  sequence  X.  by  the  method  of  cumulants,  one  is  led  to  the  study  of  the 
asymptotics  of  a  certain  type  of  deterministic  sums  associated  with  graphs.  When  the  graph  is 
a  cycle,  this  "graph-sum"  is  a  product  of  Toeplitz  matrices,  and  the  asymptotics  in  this  case 
was  obtained  by  Szego  (1958). 

Two  results  concerning  the  asymptotics  of  these  general  "graph-sums"  have  been 
obtained:  An  inequality  which  establishes  their  order  of  magnitude,  as  well  as  an  exact  limit 
theorem.  The  results  show  dependence  on  the  bond  matroid  structure  of  the  graph. 

The  main  tool  in  establishing  these  result'’  was  showing  that  a  Holder  type  inequality 
for  multiple  integrals  of  functions  applied  to  linearly  dependent  arguments  holds,  provided 
certain  rank  conditions,  known  to  physicists  as  the  "power  counting"  conditions,  are  satisfied. 

Computing  the  cumulants  by  the  diagram  formulae  and  then  applying  our  asymptotic 
results  may  reduce  the  establishing  of  rather  complicated  central  limit  theorems  to  easy  graph 
problems. 

The  method  assumes  formally  only  stationarity  of  the  underlying  sequence  Xj,  and  a 
certain  analytical  condition,  which  however  is  well  understood  only  for  Gaussian  or  linear  X-. 
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AGGREGATED  MARKOV  PROCESSES  WITH  NEGATIVE  EXPONENTIAL 

TIME  INTERVAL  OMISSION 

Frank  Ball 

University  of  Nottingham,  England 

We  consider  a  time  reversible,  continuous  time  Markov  chain  on  a  finite  state  space. 
The  state  space  is  partitioned  into  two  sets,  termed  open  and  closed,  and  it  is  only  possible  to 
observe  whether  the  process  is  in  an  open  or  a  closed  state.  Such  aggregated  Markov 
processes  have  found  considerable  application  in  the  modelling  and  analysis  of  single  channel 
records  that  occur  in  certain  neurophysiological  investigations.  The  aim  of  such  single  channel 
analysis  is  to  draw  inferences  concerning  the  structure  and  transition  rates  of  the  underlying 
Markov  process  from  the  observed  aggregated  process. 

A  further  problem  with  single  channel  analysis  is  that  short  sojourns  in  either  the  open 
or  closed  states  are  unlikely  to  be  detected,  a  phenomenon  known  as  time  interval  omission. 

We  present  a  complete  solution  to  the  problem  of  time  interval  omission,  in  the  situation  when 
the  length  of  minimal  detectable  sojourns  follows  a  negative  exponential  distribution  with 
mean  p~  .  We  provide  an  asymptotic  analysis  as  p  tends  to  infinity,  which  allows  us  to 
discuss  the  robustness  of  the  above-mentioned  structural  inferences  to  time  interval  omission. 
We  illustrate  the  theory  with  a  numerical  study  of  a  model  for  the  gating  mechanism  of  the 
locust  muscle  glutamate  receptor. 
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LOCAL  ASYMPTOTIC  MIXED  NORMALITY  OF  LOG-LIKELIHOOD 
BASED  ON  STOPPING  TIMES  AND  ITS  RATE  OF  CONVERGENCE 


A.  K.  Basu  and  Debasis  Bhattacharya 
Calcutta  University 

Let  0,  the  parameter  space,  be  an  open  subset  of  R^,  k  2  1.  For  each  0  €  0  let  the 

r.v.'s  Xm>  m  *  1,2,...  be  defined  on  the  probability  space  (£LA»Pg)  and  take  values  in  the 

space  (S,Q,  where  S  is  a  Borel  subset  of  a  Euclidean  space  and  £  is  the  o-field  of  Borel 

subsets  of  S.  It  is  assumed  that  the  joint  probability  law  of  any  finite  set  of  such  r.v.’s 

{Xm,m  2  1}  has  some  known  functional  form  except  the  unknown  parameter  0.  For  h  e  !R^ 

and  a  sequence  of  p.d.  matrices  8n  =  8n  (0g),  set  0Q  =  0  =  8q  +  5n  h,  where  0g  is  the 

true  value  of  0,  as  one  value  of  0.  For  each  n  2  1,  let  vQ  be  stopping  times  defined  on  the 

* 

process,  with  some  desirable  properties.  Let  Am(0  ,0q)  be  the  log-likelihood  ratio  of  the 
probability  measure  P  *  wxt  the  probability  measure  Pmg  ,  where  Pmg,  is  the 

restriction  of  Pq  over  =  a  (Xj^.-.-vX^.  Replacing  m  by  vn  in  Am(0  ,0g)  we 

* 

get  the  randomly  stopped  log-likelihood  ratio,  namely  Ay  (8  ,0g).  The  main  purpose  of  this 

n 

paper  is  to  show  that  under  certain  regularity  conditions  the  limiting  distribution  of 

* 

Ay  (0  ,0q)  is  locally  asymptotically  mixed  normal.  Two  examples  are  also  taken  into 
n 

account.  First  a  Berry-Essdn  Theorem  is  obtained  under  LAN  condition.  Then  a  rate  of 
convergence  is  obtained  under  LAMN  conditions. 
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AN  APPROXIMATE  MARTINGALE  CENTRAL  LIMIT  THEOREM  IN  HILBERT  SPACE 

Michael  S.  Bingham 
University  of  Hull,  U.K. 

Let  H  be  a  real  separable  Hilbert  space  with  inner  product  (•,•)•  By  an  S-operator 

we  mean  a  nonnegative  linear  operator  on  H  with  finite  trace.  Suppose  that  we  are  given  an 

adapted  triangular  array  {Snj,  <5^:  0<j<kn,n>l}  of  H-valued  random  variables  defined  on 

the  probability  space  (0,^?P).  Define  X  -:  =  S  •  —  S  .  Assume  the  following:. 

nj  nj  n,j  i 

For  every  neighbourhood  N  of  0  in  H] 

P(XRj  t  N  for  some  j  =  l,2,...Jcn)  -*  0  as  n  -»»  . 

There  is  a  neighbourhood  M  of  0  in  H  such  that  Xnj:  -  'M(VXnj  sadsfy 

kn 

l  I EKXV.y) |  ^ j  j ] |  -E-  0  as  n  -* « 
j=l 

and 

kn 

I  <Xnj,y)2  (Sy'y)  “  "  "  “ 

j=l 

for  every  yeH  where  S  is  a  random  S-operator  on  H.  Assume  also  that  there  is  a  fixed 
S-operator  T  on  H  such  that 

kn 

E <1  |E[(X;j,y)|^nj_1]|)^V(Ty^  for  all  yeH 

j=l 

and 

kn 

E[  l  *  (Ty’y)  for  a11  y  6  H 

j=l 

both  hold  for  all  sufficiently  large  n.  Finally,  assume  also  that  the  filtrations  are  nested:  i.e. 
s  ^n+lj  foraU  "J' 

Then  S  .  converges  stably  in  law  as  n  -*  »  to  the  mixture  of  Gaussian  distributions 
nKn 

with  characteristic  function:  y  -*  E[exp[-l/2(Sy,y)]],  yeH. 
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STOCHASTIC  MODELS  FOR  THE  EVOLUTION  OF  SOCIAL  BEHAVIOUR 
BASED  ON  THE  RESOURCE  DISPERSION  HYPOTHESIS 

Paul  Blackwell 

Department  of  Mathematics 
The  University  of  Nottingham 

The  resource  dispersion  hypothesis  suggests  that  heterogeneity  of  the  distribution  of 
food  may  sometimes  create  opportunities  for  animals  to  share  territories,  without  any 
cooperation  or  social  interaction.  The  evolutionary  model  presented  here  attempts  to 
determine  when  the  strategy  of  exploiting  such  an  opportunity  will  be  a  successful  one. 

A  juvenile  forming  a  group  with  its  parents  may  improve  its  chance  of  survival,  it  may 
contribute  to  its  parents'  future  reproductive  success,  and  it  may  inherit  its  parents'  territory  on 
their  death.  The  alternative  strategy  is  for  a  juvenile  to  try  to  establish  a  territory  in  its  own 
right,  and  breed  immmediately. 

The  success  of  each  strategy  is  clearly  density  dependent.  In  addition,  a  distinction 
must  be  made  between  established  adults  and  juveniles,  and  birth  and  death  rates  depend  on 
the  grouping  of  the  population.  We  model  the  system  as  a  Markov  chain,  with  the  state 
depending  not  only  on  the  number  of  users  of  each  strategy,  but  also  on  their  organization  into 
groups.  When  one  strategy  is  rare  the  system  can  be  approximated  by  a  branching  process. 
The  global  behaviour  of  the  system  can  be  investigated  using  numerical  or  simulation 
techniques. 
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LAW  OF  LARGE  NUMBERS  AND  CENTRAL  LIMIT  THEOREM  FOR  A  STOCHASTIC 
MODEL  OF  A  NONLINEAR  CHEMICAL  REACTION  WITH  DIFFUSION 


Douglas  Blount 

University  of  Utah 
Salt  Lake  City,  Utah 

A  stochastic  model  of  a  nonlinear  chemical  reaction  with  diffusion  is  constructed  by 

dividing  the  unit  interval  into  N  cells,  allowing  particles  within  each  cell  to  reproduce  as  a 

nonlinear  birth-death  process,  and  coupling  'Us  by  diffusion  with  jump  rate  proportional  to 
2 

N  .  Cell  numbers  are  divided  by  1,  a  parameter  proportional  to  the  initial  number  of  particles 
in  a  ceU.  The  rescaled  quantities  are  viewed  as  concentrations  and  the  resulting  space-time 
Markov  process  is  compared  to  the  corresponding  deterministic  model,  a  solution  to  a 
nonlinear  partial  differential  equation.  Letting  1  and  N  approach  infinity  in  an  appropriate 
manner  allows  one  to  prove  a  law  of  large  numbers  and  central  limit  theorem.  By  using  a  new 
method  the  limits  are  shown  to  hold  in  the  best  possible  state  spaces,  improving  previous 
results  of  Kotelenez. 
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A  FLUCTUATIONS  LIMIT  FOR  SCALED  AGE  DISTRIBUTIONS 
AND  WEIGHTED  SOBOLEV  SPACES 

A-  Bose 

Department  of  Mathematics  and  Statistics 
Carleton  University 
Ottawa,  Canada  K1S  5B6 

It  is  shown  that  under  the  central  limit  scaling,  the  fluctuations  of  the  space-time 
renormalized  age  distributions  of  particles  (whose  development  is  controlled  by  critical  linear 
birth  and  death  processes)  around  the  law  of  large  numbers  limit  converge  in  a  Hilbert  space 
(dual  of  a  weighted  Sobolev  space  containing  the  class  of  signed  Radon  measures  with  finite 
moment  generating  functionals)  to  a  continuous  Gaussian  process  satisfying  a  Langevin 
equation.  So  far,  the  space  of  rapidly  decreasing  functions  has  been  considered  to  be  the 
natural  state  space  for  the  kind  of  limit  theorem  considered  here.  However,  this  space  of 
rapidly  decreasing  function  is  not  suitable  in  the  present  context  and  we  are  led  to  define  a 
family  of  Sobolev  spaces  which  are  appropriate  for  our  purpose. 
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DEPENDENCIES  IN  BRANCHING  PROCESSES 


Per  Broberg 

Chalmers  Univ.  Technology 
S-41296  Gbteborg 
Sweden 

In  my  talk  I  will  survey  some  work  on  dependencies  in  branching  processes,  mainly  my 
own,  cf.  the  References. 

A  model  allowing  for  reproductive  dependencies  between  siblings  will  be  presented, 
put  in  a  wider  context  and  related  to  other  authors'  works,  including  those  of  Jagers,  Neveu  and 
Mode. 

References 

Critical  Branching  Populations  with  reproductive  Sibling  Correlations,  Stock.  Proc.  Appl.,  30 
(1988)  133-147. 

A  Note  on  the  Extinction  Probability  of  Branching  Populations,  Scand.  J.  Statist.,  14  (1987) 
125-129. 
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LIMIT  THEOREMS  VIA  LARGE  DEVIATIONS 


Wlodzimierz  Bryc 

Department  of  Mathematics 
University  of  Cincinnati 
Cincinnati,  OH  45221 


Limit  theorems,  which  can  be  obtained  from  the  Large  Deviation  Principle  will  be 
discussed.  Some  known  and  new  results  on  Large  Deviation  Proofs  of  the  Law  of  Large 
Numbers,  Central  Limit  Theorem  and  Conditional  Limit  Theorems,  each  under  appropriate 
additional  assumptions,  will  be  presented. 
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TRAPEZOIDAL  MONTE  CARLO  INTEGRATION 


* 

Elias  Masry 

Department  of  Electrical  and  Computer  Engineering 
University  of  California  at  San  Diego 
La  Jolla,  CA  92093-0407 

Stamatis  Cambanist 

Department  of  Statistics 
University  of  North  Carolina 
Chapel  Hill,  NC  27599-3260 


The  approximation  of  weighted  integrals  of  random  processes  by  the  trapezoidal  rule 
based  on  an  ordered  random  sample  is  considered.  For  processes  which  are  once  mean-square 
continuously  differentiable  and  for  weight  functions  which  are  twice  continuously 
differentiable,  it  is  shown  that  the  rate  of  convergence  of  the  mean— square  integral 

A 

approximation  error  is  precisely  n  ,  and  the  asymptotic  constant  is  also  determined. 


AMS  Subject  Classification  (1980):  65D30,  60G12,  65U05,  62G05. 

Key  Works:  Monte  Carlo  integration  of  random  processes,  trapezoidal  rule,  rate  of 
quadratic-mean  convergence. 

*Research  supported  by  the  Office  of  Naval  Rresearch  under  Contract  No. 

N00014— 84— K-0042. 

tResearch  supported  by  the  Air  Force  Office  of  Scientific  Research  Contract  No.  F49620  85C 
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MARKOV  BRANCHING  PROCESSES  WITH  INSTANTANEOUS  IMMIGRATION 


Anyue  Chen  and  Eric  Renshaw 

Department  of  Statistics 
University  of  Edinburgh,  UK 

Markov  branching  processes  with  instantaneous  immigration  (the  conditional 
expectation  of  the  sojourn  time  at  extinction  is  zero),  form  a  structure  which  represents  a 
departure  from  general  branching  processes.  In  this  paper  we  prove  that  if  the  sum  of  the 
immigration  rates  is  finite  then  no  such  structure  can  exist;  a  necessary  and  sufficient  condition 
for  existence  is  then  given  for  the  case  in  which  the  sum  is  infinite.  Study  of  the  uniqueness 
problem  shows  that  for  honest  processes  the  solution  is  unique. 
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FINITE  REVERSIBLE  NEAREST-PARTICLE  SYSTEMS 
IN  INHOMOGENEOUS  AND  RANDOM  ENVIRONMENTS 


Dayue  Chen 

University  of  California,  Los  Angeles 

In  this  paper  we  propose  and  study  finite  reversible  nearest-particle  systems  in 
inhomogeneous  and  random  environments.  Using  the  Dirichlet  principle  and  the  ergodic 
theorem  we  prove  that  a  finite  reversible  nearest-particle  system  in  random  environments 
survives  if  E  log  X  >  0  and  dies  out  if  EX  <  1.  Some  discussion  is  provided  to  show  both 
survival  and  extinction  may  happen  when  E  log  X  <  0  and  EX  >  1. 

Research  supported  in  part  by  NSF  Grant  DMS  86-01800 
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ON  SOME  INEQUALITIES  OF  LOCAL  TIME 


Chou  Chiug— Sung 

Institute  of  Mathematics 
National  Central  University,  Chung-Li, 
Taiwan,  R.O.C. 


Let  (Q,F,(Ft)t£Q,P)  be  a  probability  space  satisfying  the  usual  conditions.  Let 
L:  (Q,F)  — ♦  R+  be  a  positive  random  variable  (not  necessarily  a  (Ft)  -  stopping  time).  Let 
be  the  o-field  generated  by  Fm  and  L.  Consider 

=  {B  e  G^;  3Bt  6  Ft  such  that  Bt  n  (t< L)  =  B  n  (t<L)); 
then  (G*p  satisfies  the  following  properties: 

a)  FtcG^,t£0 

b)  (G^)  verifies  the  usual  conditions 

c)  L  is  a  (G^)  -  stopping  time. 

Let  |  Fj],  =  (Z^)  is  a  potential  of  class  (D).  The  decomposition  of 

Doob-Meyeris  Z^  =  M-A.  We  put  IT  =  inf  Z^  and  use  Z  to  denote  zK 

^  (Ks<L  s 

Our  principal  results  are  local  time  inequalities  stopping  at  any  time,  which  generalize 
results  of  Barlow— Yor. 

Theorem.  Assume  that  M  is  continuous  and  ((— —  M)^  >  is  a.s.  finite.  Then,  there  exists 
a  constant  Cp  depending  only  on  p,  1  £  p  <  such  that  for  every  continuous  (F{)  -  local 
martingale  X,  X*  the  local  time  of  X,  a  e  1R,  we  have 

1)  EpJ'P]  <;  C  E[(l+log4-)p/2  <X)P/2)1  ,  X*  =  sup  X * 

L  PI  AL  L  •>  1  aeR  1 

2)  E[XfP]  <;  C  E[(l+log— r— )p  X*P]  ,  X*  =  sup  |X  | 

L  PL  lL  L  J  1  s£t  s 

3)  E[X*P]  <;  Cp  E  [{ l+aog-4j-)p)^P] 

4)  E[(X>P/2]  <!  Cp  E  [{ l+(log-4-)P}^P] 
where  Cp  can  change  from  place  to  place. 
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UPON  A  POSSIBILITY  OF  GENERATION  OF  RANDOM  GRAPHS 

Dr.  Moise  Cocan 

Department  of  Mathematics 
University  of  Brasov 

The  present  paper  studies  a  possibility  of  generation  of  random  graphs  by  using  the 
concept  of  Markov  process. 

We  study  the  Markov  process  (^(t))t6jj  where  £(t)  =  Gt>  t  e  N,  Gt  =  (Vt,Mt),  V  *  <j>, 
|Vt|  =nt,  | Mt |  =mt,  nt  +  mt  =  t,  te  IN,  and  Gx  =  (VpCtO,  jVJ  =  1. 

Let's  presume  that  Gt  is  obtained  from  Gt_j  by  adding  a  new  node  or  a  new  edge. 

It  is  clear  that  if  the  graph  G£_^  is  complete,  Gt  is  obtained  by  adding  a  new  node.  The 
distribution  of  probability  of  the  random  variable  Gt  (Gt  is  a  variable  with  graph  values)  is 

(Vt_i  u  (v),  Mj_j)  Vt-l’Mt-l  u  {m,) 

p  q 

where  p,q  >  0,  p  +  q  =  1,  v  i  V{_p  m  €  M£_p 

Different  feature  of  graphs  Gt(t€lN)  are  studied  and  a  certain  algorithm  is  given  for  the 
generation  of  sample  subgraphs  of  the  graph  Gt(teW). 

For  the  generation  of  graphs  Gt(telN),  an  algorithm  is  given  whic..  is  able  to  select 
those  that  have  certain  features. 

Considering  graph  Gt(te!N),  t-fixed,  as  a  structured  population  and  the  sample 
subgraphs  as  samples  from  this  population,  a  study  of  possibilities  of  estimation  of  parameters 
that  show  the  character  of  the  considered  graph  is  realised. 
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HITTING  PROBABILITIES  FOR  SPECTRALLY  POSITIVE  L&VY  PROCESSES 

R.  A.  Doney 
York  University 

Let  (Y(t),t  SO}  be  a  spectrally  positive  Livy  process,  that  is  a  process  with  stationary 
and  independent  increments  which  has  no  negative  jumps,  and  let  x(x)  denote  the  hitting 
probability  of  the  singleton  {x},  viz  x(x)  =  Pr{Y(t)  =  x  for  some  0  <  t  <  «}.  We  give  a 
formula  in  terms  of  the  Livy  exponent  of  Y  for  the  Laplace  transform  of  x(x),  and  use  it  to 
determine  explicitly  the  value  of  x(0)  and  i  im  x(x).  In  the  case  that  Y(t)  =  U(t)  - 1,  where 

X->» 

U  is  a  subordinator,  or  non-decreasing  Levy  process,  the  first  result  confirms  a  recent 
conjecture  of  Mallows  and  Nair  and  the  second  extends  a  result  of  Nair,  Shepp  and  Klass  for 
the  case  that  U  is  a  Poisson  process.  We  also  give  formulae  for  the  Laplace  transform  of 
ECe-®^^)  and  E(e~®D^),  where  T(x)  is  the  time  at  which  the  first  visit  to  x  takes 
place,  and  D(x)  is  the  time  between  the  first  passage  over  x  and  T(x). 
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DUALITY  PRINCIPLE  IN  DOUBLY  CONTROLLED  PROCESSES  OF  SERVICING 
MACHINES  WITH  RESERVE  REPLACEMENT 

Jewgeni  Dshalalow 

Department  of  Applied  Mathematics 
Florida  Institute  of  Technology,  USA 

The  classical  system  of  servicing  m  +  1  unreliable  machines  by  a  single  repairman 
was  studied  by  Takdcs  (1962).  In  the  present  modification  of  this  model,  the  author  allows  w 
additional  machines  to  be  placed  on  reserve  substituting  any  failed  one.  In  addition,  a 
comprehensive  control  over  the  failure  rates  and  repair  times  maintains  the  whole  system.  The 
author  introduces  two  models,  with  and  without  idle  periods  of  the  repairman,  and  he  defines 
and  applies  a  duality  principle  to  find  an  explicit  relation  between  the  servicing  processes  in 
both  models.  The  author  further  develops  semi-regenerative  techniques  (invented  earlier, 
1985,  1989a)  and  applies  them  to  the  servicing  process  in  the  model  without  idle  periods 
which  seems  to  be  simpler  analytically.  Then,  the  duality  relation  allows  to  extend  the  results 
obtained  for  a  more  attractive  system  with  idle  period 
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THE  DISTRIBUTION  OF  THE  FREQUENCIES  OF  AGE-ORDERED  ALLELES 

IN  A  DIFFUSION  MODEL 


S.  N.  Ethier 
University  of  Utah 

We  prove  that  the  frequencies  of  age-ordered  alleles  in  the  stationary 
infinitely— many— neutral— alleles  diffusion  model  are  distributed  as 

x1,(i-x1)x2,(i-x1xi-x2)x3,..., 

where  X^^^,...  are  independent  beta(  1,9)  random  variables,  9  being  twice  the  mutation 
intensity;  that  is,  the  frequencies  of  age-ordered  alleles  have  the  so-called 
Griffiths-Engen-McCloskey,  or  GEM,  distribution.  In  fact,  two  proofs  are  given,  the  first 
involving  reversibility  and  the  size-biased  Poisson-Dirichlet  distribution,  and  the  second 
relying  on  a  result  of  Donnelly  and  Tavard  relating  their  age-ordered  sampling  formula  to  the 
GEM  distribution. 
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CYCLIC  CELLULAR  AUTOMATA  AND  RELATED  PROCESSES 


Robert  Ksch 

University  of  North  Carolina  at  Charlotte 

A  cyclic  cellular  automaton  is  a  descrete— time  process  defined  on  the  state  space 

{0,1,...,N— 1}^.  Initially,  each  lattice  site  of  2^  is  randomly  assigned  one  of  the  values  from 
{0,1,. ...N-l},  which  are  called  colors.  The  N  colors  are  arranged  in  a  cyclic  heirarchy;  we 
say  that  color  i  can  eat  color  j  if  i— j=l  (mod  N).  Given  an  initial  configuration,  the 
evolution  proceeds  deterministically.  At  each  time  step,  each  site  x  looks  at  the  color  of  each 
of  its  neighbors.  If  the  color  of  any  neighbor  can  eat  the  color  of  x,  then  it  does  so,  meaning 
that  x  changes  its  color  to  the  eating  color.  All  sites  update  synchronously  in  this  manner  at 
each  time  step. 

Cyclic  cellular  automata  have  properties  akin  to  cyclic  particle  systems,  introduced  by 

Bramson  and  Griffeath,  which  are  continuous  time  Markov  processes  on  {0,1,...,N-1}  with 
a  stochastic  evolution  rather  than  a  deterministic  one.  Some  empirical  observations  of  the  two 
cyclic  systems  in  two  dimensions  will  be  presented  as  qualitative  comparisons.  Then,  the 
focus  will  shift  to  the  one-dimensional  cyclic  cellular  automaton  to  present  some  results 
concerning  clustering.  Finally,  we  shall  discuss  the  relation  of  this  clustering  result  to  other 
theorems  about  one-dimensional  systems  of  particles  with  various  interaction  mechanisms. 
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REGRESSION  MODELS  WITH  INCREASING  NUMBER  OF  UNKNOWN  PARAMETERS 


A.  G.  Gadiiev 


Institute  of  Cybernetics 
Azerbaijanian  Academy  of  Sciences 
Agaev  str.,  9,  Baku,  370105,  USSR 


We  consider  the  regression  model 


=  Tl(xi,0)  +  i  = 

T 

where  ti(Xj,0)  is  a  some  nonlinear  function,  0  =  (0j,...,0  )  is  the  vector  of  unknown 


parameters.  Assume 


E^  =  0,  E£?  =  o?  -  are  unknown  and  °i  So0 


The  number  m  of  unknown  parameters  depends  on  N  and 


m/N  -*  0,  N  -* « 


Let  0  denote  a  least  squared  estimator 


a  dnuJ) 


A  A  AAA 


F-(f*j).y =f-0.Iw(0)  =  (V 


A-.  =  5.- 

y  y 


dqCxj,©)  3ti(x.,0) 


ATA 

,F 1  Fx-1 


0  <  £  ...  £  are  eigenvalues  of  the  matrix  \  =  NE(0-0)(0-0)T  is  the 

covariated  matrix  of  the  vector  VN(0-0),  are  the  elements  of  the  matrix  CN  , 

ak£  =  (y-y)Tljc,£(9)(y-y)»  A  =  (akp,  k,£  =  I7m 

AJA  A'J'A 

and  -  is  the  element  (k,£)  of  the  matrix  CN  =  Let 

m/N’X^  ->  0  for  N-*«  (3) 

m- ,/m/N' Aj  ->  0  for  N  -♦  »  .  (4) 
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Theorem.  If  conditions  (1)— (3)  are  true  then 

(8-0)  ^  0,  N  -  .  (*) 

If  condition  (4)  is  also  true  then 

E(Ck£-ck£)  ■*  °*  (^k£_Ck£)  ^  0  for  N  -  -  .  (**) 

m 

If  m  <  c  <  <*>  and  ti(x.,0)  =  ^  0j  -^(x.),  i.e.,  we  have  linear  regression,  then 

j  =  l 

Corollary.  If  lAj  -*  0  as  N  ->  ~  then  (*)  and  (**)  hr 

A 

The  matrix  CN  can  be  used  for  the  construction  of  a  confidence  band  for  the  function 
ti(x,0).  For  example  if  £■  are  Gaussian  then  the  quadratic  form 

A  ta— 1  A 

(eVcj/^-e) 

has  chi-squared  distribution  and  it  is  easy  to  construct  a  confidence  band  to  the  function 
ti(x,0),  for  example  according  to  the  approach  suggested  [1]. 
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UNIQUENESS  OF  VISCOSITY  SOLUTIONS  (V.S.)  OF  HAMILTON— JACOBMJELLMAN 
(HJB)  EQUATIONS  FOR  CONTROLLED  DIFFUSIONS  ON  FINITE-DIMENSIONAL 
RIEMANNIAN  MANIFOLDS  WITH  BOUNDARY 

Giuliano  Galeazzi 

via  Gaetano  Ciarrocchi,  15-00151, 

Rome,  Italy 

The  notion  of  V.S.  can  be  extended  (talk  at  17th  Conference)  to  HJB  Equations  on 
Riemannian  Manifolds  M's: 

G(p,V)  =  inf  {Aa  ( p,aQ) V(p)+f(p,aQ) )  =  0  .  (C  =  "Control  actions" 

V(p)  =  0  on  9M  .  Aa  =  "Different.  Gen.") 

in  a  suitable  way.  Uniqueness  is  now  proven  by  semi-group  arguments.  A  class  of  second 
order  elliptic  operators  "A'"s  has  the  property  (R.  Azencott)  that  (Pp,  the  minimal  strong 
Markov  semigroups  associated  with  the  A's,  leave  Cq(M)  (C^(M)  &  vanishing  at  ») 
invariant  and  turns  Borel  &  Bounded  g's  into  Ptg  e  (^(M).  As  M  is  compact  and  V.S.  are 
C°(M),  this  holds  for  differences  of  any  "two"  V.S.:  P^Vj-V^  €  C^(M).  This  extended 
notion  of  V.S.  implies  the  existence-on  geodesic  balls  Bp  centered  at  point  p  6  M  —  of  two 
subsets  D+V(p)  and  D-V(p)  of  C^Bp)  that  have,  respectively:  G(p,H)  ^  0,  for 
H  6  D+V(p)  a~d  G(p,H)  £  0  for  HeD~V(p).  Now  if  the  Aa  (elliptic)  obey  the  Positive 
Maximum  Principle  (the  A's  automatically  do),  have  A(p)g(p)  £  Aa(p)g(p)  for  g  e  C^(M) 
and,  finally,  a  local  condition  on  A's  minimal  semigroup  (P£)  holds,  then  for  any  "two"  V.S.: 
A(Pt(Vj-V2))  £  Aa(HI-H2)  as  Hj  e  D+V^p),  Rj  e  D~V2(p).  Under  the  existence  of  a 
dense  subset  S  of  M  at  whose  points  those  spaces  "D"  are  nonvoid  for  all  V.S.  such  an 
inequality  gives  A^V^-V^))  £  0  on  M,  to  the  effect  that  the  Classical  Maximum  Principle 
yields  P^Vj-V^)  £  0  (t  sufficiently  small).  (P^'s  being  strong  entails  Vj-V2  £  0.  Finally, 
exchanging  Vj's  and  V2's  roles  yields  the  converse  inequality. 
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A  GEOMETRIC  RATE  OF  CONVERGENCE  FOR  THE  MAXWELL  GAS 
WITH  AN  ANGULAR  CUT-OFF 


Brigitte  Chauvin 
University  Paris  VI 
Laboratoire  de  Probability 
4,  Place  Jussieu,  tour  56 
75230  Paris  Cedex  05 

Gaston  Giroux1 

Ddpartement  de  Mathdmatiques 
University  de  Sherbrooke 
Sherbrooke,  Quebec 
J1K2R1 


We  consider  the  measure-  alued  time-dependent  solution  to  the  Boltzmann  equation 
for  a  Maxwell  gas  with  an  angular  cut-off.  For  this  solution,  we  show  that,  for  any  initial  law 
having  all  its  moments  finite  up  to  some  integer,  convergence  to  the  equilibrium  along  good 
test  functions  occurs  at  a  geometric  rate. 

Our  proof  is  based  on  a  representation  of  an  underlying  process  involving  an 
exponential  age-dependent  binary  branching  process.  Some  careful  computation  of 
conditional  expectation  is  also  needed. 


This  gives  an  answer  to  an  old  problem;  see  for  instance: 


McKean,  H.  P.  Speed  of  approach  to  equilibrium  for  Kac's  caricature  of  a  Maxwellian  gas. 
Arch.  Rat.  Mech.  Anal.  21  (1966)  343-367. 

Arkeryd,  L.  Stability  in  L1  for  the  space  homogeneous  Boltzmann  equation.  Technical 
Report.  University  of  Gdteborg,  Sweden  (1988). 

Research  supported  in  part  by  NSERC,  Canada,  under  Grant  No.  A-5365. 
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NONLINEAR  LIMIT  FOR  A  SYSTEM  OF  DIFFUSING  PARTICLES 
WHICH  ALTERNATE  BETWEEN  TWO  STATES 

Carl  Graham 

CMAP,  Ecolc  Polytechnique 
91128  Palaiseau  FRANCE  (URA  CNRS  756) 

We  consider  a  system  of  interacting  particles,  each  evolving  alternately  under  two 
different  diffusion  dynamics  in  1R^.  For  example,  these  particles  could  be  alternately  in  a 
gaseous  and  in  a  liquid  phase,  or  adsorbed  and  desorbed  by  a  medium.  A  special  case  is 
McKean's  2-speed  model  for  the  Boltzmann  equation. 

We  study  the  corresponding  nonlinear  McKean- Vlasov  diffusion,  where  a  particle  sees 
the  medium  as  a  function  of  its  law.  This  can  be  considered  as  the  Statistical  Mechanics 
weak-interaction  limit  for  the  system,  through  a  propagation  of  chaos  result. 

This  is  a  Probabilistic  model  for  a  Chromatography  tube,  where  molecules  in  a  gaseous 
phase  pushed  by  a  flow  of  neutral  gas  are  adsorbed  and  desorbed  by  a  liquid  phase  in  the  tube. 
The  time  taken  by  the  molecules  to  get  through  the  tube  is  measured,  and  the  device  can  be 
used  to  identify  a  substance,  separate  the  components  of  a  mixture,  and  so  on. 

Nonlinear  problems  may  often  be  dealt  with  by  contraction  methods.  This  is 
straightforward  in  a  stochastic  differential  equation  setting,  or  for  a  pure-jump  process  whose 
generator  integral  form  allows  to  work  directly  on  the  martingale  problem.  The  problems  here 
arise  from  the  nonexistence  of  a  good  strong  pathwise  representation  and  from  the 
integro-differential  generator. 

Our  idea  here  is  to  construct  processes  having  the  appropriate  laws,  choosing  the  joint 
law  so  that  the  sample  paths  are  as  close  as  possible,  and  then  use  stochastic  calculus.  This  is 
automatically  done  for  strong  stochastic  differential  equations  by  taking  the  same  driving 
terms.  Here  we  use  time-change  first,  in  order  to  get  to  a  situation  where  strong  pathwise 
contraction  techniques  apply.  This  is  a  rather  new  technique,  developed  in  a  previous  joint 
paper  with  Professor  M6tivier.  Thus,  we  obtain  existence  and  uniqueness  for  the  nonlinear 
process.  Then  we  prove  propagation  of  chaos. 
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MAX  INFINITELY  DIVISIBLE  AND  MAX  STABLE  SAMPLE 
CONTINUOUS  PROCESSES 


Evarist  Gind 

CUNY (The  College  of  Staten  Island)  and  Texas  A&M  University 

Marjorie  G.  Hahn 
Tufts  University 

Pirooz  Vatan 
Strategic  Automation  Inc. 


Conditions  for  a  process  5  on  a  compact  metric  space  S  to  be  simultaneously  max 
infinitely  divisible  and  sample  continuous  are  obtained.  Although  they  fall  short  of  a  complete 
characterization  of  such  processes,  these  conditions  yield  complete  descriptions  of  the  sample 
continuous  non-degenerate  max  stable  processes  on  S  and  of  the  infinitely  divisible 
non— void  random  compact  subsets  of  a  Banach  space  under  the  operation  of  convex  hull  of 
unions. 
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CORRELATION  THEORY  OF  ALMOST  PERIODICALLY  CORRELATED  PROCESSES 


Hairy  L.  Hurd 

Center  for  Stochastic  Processes 
Department  of  Statistics 
University  of  North  Carolina 
Chapel  Hill,  N.C.  27599-3260 


A  second  order  stochastic  process  X  is  called  almost  periodically  correlated  (PC)  in 
the  sense  of  Gladyshev  if  its  mean  function  m(t)  and  covariance  R(t+x,t)  are  almost  periodic 
functions  of  t  for  every  x.  We  show  that  the  mean  uniformly  almost  periodic  processes 
discussed  by  Kawata  are  also  almost  PC  in  the  sense  of  Gladyshev.  If  X  is  almost  PC, 
then  for  each  fixed  x  the  function  R(t-rX,t)  has  the  Fourier  series  representation 


R(t+x,t)  -  V  a(X,x)exp(iXt) ,  and 
A+c 

i  r 

a(X,x)  =  1  im  ££  R(t+x,t)exp(-iXt)dt 

A”*00  ^  A 

-A+c 


(0 


(u) 


exists  for  every  X  and  x,  independently  of  the  constant  c.  Assuming  only  that  a(X,x) 
exists  in  this  sense  for  ewry  X  and  X,  we  show  a(X,x)  is  a  Fourier  transform 


a(X,x)  = 


!k 


exp(iyr)r^(dy) 


(iii) 


if  and  only  if  a(0,x)  is  continuous  at  x  =  0;  under  this  same  condition,  the  set  A  =  (X: 
a(X,i)  *■  0  for  some  x}  is  countable.  We  show  that  a  strongly  harmonizable  process  is  almost 
PC  if  and  only  if  its  spectral  support  is  contained  in  {(YpY^'  Y2  =  Y^  -  ^  «  A},  a 

countable  set  of  straight  lines  parallel  to  and  including  the  main  diagonal;  further,  one  may 
identify  the  spectral  measure  on  the  k^1  line  with  the  measure  r,  appearing  in  (iii).  Finally 


we  observe  that  almost  PC  processes  are  asymptotically  stationary  and  give  conditions  under 
which  strongly  harmonizable  almost  PC  processes  may  be  made  stationary  by  an  independent 
random  time  shift. 
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ON  THE  EXCHANGEABILITY  OF  DIFFERENTIATION  AND  STOCHASTIC 
INTEGRATION  WITH  APPLICATIONS  TO  STOCHASTIC  PDE'S 


Jing— Lin  Jiang 

Department  of  Mathematics 
Wayne  State  University 
Detroit,  MI  48202 


In  this  paper,  we  establish  a  result  on  the  exchangeability  of  differentation  and 
stochastic  integration,  namely, 


T 


d  f 
3x 

0 


f(s,x,co)dWs 


T 

’  3*f<s-*’“>dws 
0 


under  some  condition  on  f.  A  similar  proposition  is  also  true  fcr  lgher-order  derivatives. 
The  results  will  be  applied  to  the  following  two  types  of  stochastic  partial  differential 
equations  to  yield  classic  solutions 


au 

3T 

=  LLU+  [f(t,x,co),Wt] 

in  D 

U(0 

,  X  ,  (0  )  =  <J)(x,0>) 

in  D 

(1) 

U(t 

,x,co)  =0 

on  3D 

au 

cTT 

=  LU  +  [hU  +  C-VU,  Wt] 

in 

(2) 

U(0 

i  X  ,  CO)  =  (j)  (  X  ,  03) 

in  1R11 

where  L  is  a  second  order  strong  elliptic  operator,  (Q,7,P,7t,Wt)  is  an  IR^— valued  standard 
Weiner  process,  D  is  a  bounded  domain  in  Rn,  f:  [0,TJ  x  d  x  Q  ->  IR^,  h:  [0,T]  x  IRn  -» IR^, 
c:  [0,T]  x  lRn  -  IR^0,  [•  ,•]  denotes  the  scalar  product  in  IR^.  An  example  of  (1)  is  the  heat 
equation  (L  =  A,  the  Laplacian)  with  the  heat  source  f(t,x,co)Wt  of  white-noice  type.  (2) 
arises  from  the  nonlinear  filtering  problem  of  diffusion  processes. 
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RANDOM  TIME  CHANGE  AND  INTEGRAL  REPRESENTATIONS 
FOR  MARKED  STOPPING  TD<ES 

Olav  Kallenberg 
Auburn  University 

The  distributions  of  all  random  triples  (x,u,T|),  such  that  (t,u)  is  a  marked  stopping 
time  with  compensator  tj,  form  a  convex  set  M,  whose  extreme  points  are  the  distributions 
obtained  when  the  filtration  is  the  one  induced  by  (x,u).  Moreover,  every  element  of  M  has 
a  unique  integral  representation  as  an  average  over  extreme  points. 

Using  the  Doldans  exponential  process  associated  with  q,  one  may  define  a  certain 
discounted  version  C,  of  the  compensator,  and  it  is  shown  that  if  T  is  a  predictable  mapping 
into  [0,1]  such  that  £T— ^  is  bounded  by  Lebesgue  measure,  then  T  is  uniformly 
distributed  on  [0,1].  The  proof  of  this  uses  the  peculiar  fact  that  EV  =0  for  any 
predictable  process  V  satisfying  a  certain  moment  condition. 

By  iterating  the  mentioned  time  change  result,  one  obtains  a  reduction  of  an  arbitrary 
simple  point  process  £  to  Poisson,  through  a  predictable  mapping  based  on  the  compensator 
Tj,  together  with  suitable  independent  randomizations,  needed  to  resolve  the  possible 
discontinuities  of  r\  at  the  atoms  of  Mention  will  also  be  made  of  a  very  general  version 
of  Knight's  and  Meyer's  time  change  reduction  of  orthogonal  continuous  local  martingales  or 
quasi-leftcontinuous  point  processes  to  independent  Brownian  motions  or  Poisson  processes, 
respectively. 
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SEMIMARTINGALES  FORCED  ONTO  A  MANIFOLD  BY  A  LARGE  DRIFT 


G.  S.  Katzcnbcrger 
University  of  Wisconsin 

We  consider  a  sequence  of  R^— valued  semimartingales 

t 

Xn(t)  "  Zn(t)  +  /  F(Xn(s-))dAn(s), 

0 

where  {Zn)  a  "well-behaved"  sequence  of  semimartingales,  F  is  a  vector  field  whose 
deterministic  flow  has  an  asymptotically  stable  mar  >ld  of  fixed  points  T  and  An  is  a 
nondecreasing  process  that  asymptotically  puts  infinite  mass  on  every  interval.  This  situation 
occurs  frequently  in  deriving  diffusion  approximations  for  certain  genetics  problems,  branching 
processes  and  stochastic  approximation  problems.  Intuitively,  if  XR  starts  close  to  T  then 
the  drift  term  forces  Xn  to  stay  close  to  T  and  the  limiting  process  (if  there  is  one)  must 
aciuaily  stay  ou  T,  hence  is  typically  lower  dimensional  than  the  original  Xn's.  We  give 
conditions  under  which  {XnJ  is  relatively  compact  in  the  Skorohod  topology  and  obtain  a 
stochastic  differential  equation  for  the  limiting  process(es) 
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A  VISCOSITY  ALGORITHM  FOR  THE  SCALAR  CONSERVATION  LAW 

Andrzej  Korzeniowslri 

Department  of  Mathematics 
University  of  Texas  at  Arlington 
Arlington,  TX  76019 

Consider  a  parabolic  equation 

2 

3u  ,  ,3  u  3u 

<Jt“  =  +  B(uo)3x_  •  uo(0’x)  =  g(x)  (1) 

which  reduces  to  a  nonlinear  hyperbolic  equation  known  as  the  conservation  law,  when  the 
viscosity  a  is  0.  Unlike  the  standard  numerical  schemes,  which  approximate  differential 
operators  by  difference  operators,  we  devise  an  algorithm  for  solving  (1)  that  utilizes  the 
Feynmann-Kac  representation  of  the  solution  u^.  It  is  well  known  that  as  a  ->  0  then 
ua-u  which  is  a  solution  to 

^=B(u)^,u(0,x)  =  0.  (2) 

We  compare  our  method  with  the  Glimm’s  scheme,  developed  for  (2),  by  means  of  numerical 
calculations  for  a  concrete  example.  Our  scheme  is  deterministic  as  opposed  to  Glimm’s 
scheme,  which  requires  a  simulation  of  uniform  random  variables,  and  as  such  suffers  no 
random  number  generation  error. 
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THE  HEAVY  TRAFFIC  APPROXIMATION 
FOR  CLOSED  EXPONENTIAL  QUEUEING  NETWORK 


E.  V.  Krichagina 


Institute  of  Control  Sciences 
Profsoyuznaya  65 
117806  Moscow  GSP-7 
USSR 


This  report  presents  limit  theorems  for  queue  length  processes  associated  with  a  general 
closed  queueing  network.  The  first  theorem  states  that  normalized  sequences  of  queue  length 
processes  converge  in  probability  x  deterministic  function  as  the  number  of  customers  in  the 
network  goes  to  infinity.  The  lit  /unction  describes  a  unique  solution  of  a  certain 
differential  inclusion.  The  secc  limit  theorem  establishes  weak  convergence  of  normalized 
deviations  of  the  queue  processes  from  a  deterministic  limit  to  a  diffusion  process  with  values 
in  some  convex  region  and  instantaneous  reflections  on  the  boundaries. 
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LOCAL  TIMES  FOR  SUFERDEFFUSIONS 


Steve  Krone 

University  of  Massachusetts-Amherst 

We  study  local  times  for  a  class  of  measure-valued  Markov  processes  known  as 
superdiffusions.  These  processes  arise  as  high-density  limits  of  systems  of  independent 
branching  IR^-valued  diffusions.  As  for  some  ordinary  Markov  processes,  superdiffusions 
(when  d  =  1,2,3)  have  local  times.  We  will  present  some  regularity  results  for  these  local 
times  (e.g.,  joint  continuity  and  Holder  condition)  and  implications  about  the  behaviour  of  the 
superdiffusion  (e.g.,  information  about  the  Hausdorff  dimension  of  the  "level  sets"). 
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MEASURING  THE  LONG-TERM  DEPENDENCE  OF  SOME 
STATIONARY  STABLE  PROCESSES 


Joshua  Levy 

University  of  Lowell 
One  University  Ave. 

Lowell,  MA  01854 

To  estimate  the  long-term  dependence  of  a  stationary  stable  process,  Y  =  (Y(t)}. 

<  t  <  «,  we  can  evaluate  the  function  r:  R  -»  C  given  by 

r(t):  =  rY(0lt02; t)  =  EexpU^YCt)  +  ©2Y(0))}  -  (Eexp^YtolXEexp^YfO)}) 
as  1 1 1  -  «  for  all  non-zero  real  numbers  0j  and  ©2-  We  will  study  the  case  when  Y  is 
the  increment  process  of  either  linear  fractional  Ldvv  motion  (I  FLM)  or  log-fractional  Ldvv 
motion.  LFLM  is  an  a-stable  H-self-similar  process  defined  for  0  <  a  £  2,  0  <  H  <  1, 
and  H  *  L-  while  log-fractional  Ldvy  motion,  which  is  defined  only  when  1  <  a  <  2,  is 
a-stable  and  ^  -self-similax 

One  of  the  results  for  LFLM  states  that  r(t)  ~  Altl^-**,  A  =  A(a,H,0p©2),  if 
1  <  a  <  2  and  1  —  <  H  <  1.  The  current  research  extends  earlier  work  of 

Austrauskas  (1984,  1987)  who  studied  the  process  fractional  Lgvv  motion,  which  is  one  of  the 
symmetric  LFLM’s.  The  constant  A  is  important  because  there  are  a-stable  H-self-similar 
processes,  Z  =  {Za  H(t)},  with  1  <  a  <  2  and  1  -  <  H  <  1,  that  have  stationary 

increments  and  that  are  generated  by  renewal  sequences  for  which  there  exist  0j  =  0j(a,H) 
and  ©2  =  ©2(a,H)  such  that  r^(t)  -  A'  1 1 1  a,  but  A'  *  A.  As  a  consequence,  Z  cannot 
be  any  LFLM. 
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BROWNIAN  MOTION  ON  NESTED  FRACTALS 


Tom  Lindstrgm 

Department  of  Mathematics 
University  of  Oslo,  Norway 

The  purpose  of  this  paper  is  to  construct  Brownian  motion  on  a  reasonably  general 
class  of  self-similar  fractals.  To  this  end,  I  introduce  an  axiomatically  defined  class  of  "nested 
fractals",  which  satisfy  certain  symmetry  and  connectivity  conditions,  and  which  also  are  (in 
the  physicists'  terminology)  finitely  ramified.  On  each  one  of  these  nested  fractals,  a  Brownian 
motion  is  constructed  and  shown  to  be  a  strong  Markov  process  with  continuous  paths.  If  the 
Laplacian  A  on  the  fractal  is  defined  as  the  infinitesimal  generator  of  the  Brownian  motion, 
and  n(a)  denotes  the  number  of  eigenvalues  of  -A  less  than  a,  I  prove  that 

n(a)-od  l°ev/1°eX  as  a--, 

where  d  is  the  Hausdorff  dimension  of  the  fractal,  and  v  and  X  are  two  parameters 
describing  its  self-similarity  structure.  In  general,  log  v/log  X  *  1/2  and  hence  Weyl's 
conjecture  can  only  hold  for  fractals  in  a  modified  form. 


AMS  1980  subject  classifications:  Primary  60J60,  60J65,  secondary  03H05,  35P20 

Key  words  and  phrases:  Brownian  motion  on  fractals,  diffusions,  self— similar  fractals,  fractal 
Laplacians,  as?  .nptotic  distribution  of  eigenvalues,  nonstandard  analysis. 
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USING  RENEWAL  THEORY  TO  INVESTIGATE  SELF-AVOIDING  WALKS 


Neal  Madras 


Department  of  Mathematics 
York  University 
4700  Keele  Street 
Downsview,  Ontario  M3J  1P3 
Canada 


It  is  possible  to  decompose  some  self-avoiding  walks  in  a  manner  which  translates 
questions  about  enumeration  into  questions  about  probabilities  of  renewal  processes.  (See,  for 
example,  H.  Kesr  J.  Math.  Phys.  4.  960  (1963),  or  J.  T.  Chayes  and  L.  Chayes,  Commun. 
Math.  Phys.  1Q5.  (1986)).  We  review  this  technique  and  show  how  probabilistic  reasoning 

can  be  used  to  derive  a  rigorous  (but  probably  nonoptimal)  upper  bound  on  the  critical 
exponent  for  the  number  of  self-avoiding  polygons  (also  called  self-avoiding  loops). 
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EXPONENTIAL  STABILITY  FOR  STOCHASTIC 
DIFFERENTIAL  EQUATIONS  WITH  RESPECT  TO  SEMIIMARTINGALES 


Xuerong  Mao 

Mathematics  Institute 
University  of  Warwick 
Coventry  CV4  7AL,  England 


Let  be  a  complete  probability  space  with  the  filtration  { 7} 

satisfying  the  usual  conditions,  i.e.,  which  is  right  continuous  and  contains  all  P-null  sets. 
Consider  the  following  stochastic  differential  equation  with  respect  to  semimartingales 

dX(t)  -  AX(t)dp(t)  +  G(X(t),t)dM(t),  t  S  tQ£  0)  (1) 

X(t(P  =  x0 

T 

where  X  =  (Xp...,Xn)  ,  A  is  a  constant  n  x  n  matrix,  p  is  a  continuous  {^-adapted 

T 

nondecreasing  process,  M  =  Mm)  is  a  m-dimensional  continuous  local  martingale 

with  M(0)  =  0,  xn  is  an  n-dimensional  7.  -measurable  vector,  and 
U  tQ 

G:  Rn  x  IR+  x  O.  -* R11*111  is  sufficiently  smooth  for  the  existence  and  uniqueness  of  the 
solution.  Denote  by  X(t,tQ,XQ)  the  solution  of  Eq.  (1). 

Eq.  (1)  might  be  regarded  as  a  stochastic  perturbed  system  of 

dX(t)  =  AX(t)dp(t) .  (2) 


Suppose  Eq.  (2)  is  exponentially  stable  almost  surely.  Under  what  conditions  is  Eq.  (1)  still 
exponentially  stable  almost  surely? 

Throughout  this  paper  we  will  assume  there  exist  [7^] -adapted  processes  K^f), 
i,j  =  l,2,...,m,  such  that 


<Mi,Mj)(t)=  f  K..(s)dp(s),  t  >  0 

D 


(3) 


Define 


K  •'=  <Kij  w 


In  this  paper  we  prove  the  following  main  result. 
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Theorem.  Let  XQ  be  the  maximum  of  the  real  parts  of  all  eigenvalues  of  -A.  Suppose  there 
exist  positive  constants  <Xp a2,  nonnegative  constants  Ppp2»Pi$2’  and  polynomials  p^(t), 
p2(t)  such  that 

o  (~^<xi^n+P  1  ^ 

||G(x,r)r£pA(t)e  zu  a.s.  (4) 

Po* 

!|K(t)|j  <;  p2(t)e  z  a.s.  (5) 

otjt  -  P1  2  ji(t)  5  o^t  +  (32  a.s.  (6) 


Furthermore,  assume 


limsup  lloglle^ll^  £  — X 
t-*»  1 


where  X  is  a  positive  constant.  Then 


1  2 

limsup  -  log||X(t,t0,x0) ||  S-ajX  +  ^PjTpj)  a.s. 
t  -*»  1 


for  all  tQ  >  0  and  -measurable  vector  Xq. 
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STEADY  STATE  ANALYSIS  FOR  A  SERIES 
QUEUEING  NETWORK  WITH  BLOCKING 

Alain  Jean-Marie 

INRIA  Sophia-Antipolis,  Valbonne  06565,  FRANCE 
William  A.  Massey 

AT&T  Bell  Laboratories,  Murray  Hill,  NJ  07974,  USA 

Finding  the  limiting  equilibrium  distribution  for  the  finite  state  Markov  chain 
associated  with  a  queueing  network  that  has  finite  buffers  is  a  notoriously  difficult  problem. 

In  this  paper,  we  obtain  an  exact  solution  and  new  results  for  the  equilibrium  behavior  of  a 
class  of  series  queueing  networks  with  blocking.  Each  node  has  capacity  one  or  buffer  size 
zero.  These  systems  arise  in  the  modeling  of  data  communication  links  or  manufacturing 
assembly  lines.  A  previous  paper  by  the  second  author,  used  combinatorial  methods  to 
analyze  such  a  network  with  the  single  Poisson  input  rate  equal  to  the  processing  rate  of  all  the 
N  exponentially  distributed  service  times.  We  generalize  this  model,  by  allowing  the  Poisson 
input  rate  X  to  differ  from  the  processing  rates  for  each  node,  all  of  which  are  normalized  to 
equal  one.  Aggregate  mean  value  quantities  for  this  system,  such  as  the  throughput  and  the 
total  number  in  the  system  were  solved  by  Bocarov.  Liggett  also  computed  the  throughput 
rate,  viewing  this  model  as  a  truncated  version  of  an  asymmetric  exclusion  process.  We 
extend  these  results  to  give  a  closed  form  solution  for  the  complete  joint  queue  length 
distribution,  as  well  as  the  marginal  distribution  for  any  subset  of  the  nodes. 

These  exact  formulas  are  achieved  by  obtaining  a  recursion  formula  for  the 
probabilities  when  suitably  normalized.  This  leads  to  a  simple  expression  for  the  grand 
partition  function  of  the  family  of  networks  for  fixed  X,  indexed  by  N.  We  also  use  this 
generating  function  to  develop  simple  estimates  of  these  marginal  probabilities.  Specifically, 
we  obtain  for  large  N,  simple  asymptotic  formulas  for  the  marginal  probabilities  of  a  given 
node  at  the  beginning,  middle,  and  end  of  the  network.  Finally,  we  obtain  an  exact  formula 
having  the  flavor  of  Little's  law,  that  relates  the  throughput  of  the  system  to  the  average 
number  of  contiguous  busy  nodes  ning  from  the  first  node  of  the  network. 
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U-STATISTICS  AND  THE  DOUBLE  STABLE  INTEGRAL 


loop  Mijnheer 

Department  of  Mathematics 
University  of  Leiden 
The  Netherlands 


Multiple  Wiener  integrals  appear  in  the  theory  of  U-statistics.  See[l].  Several  results 
are  known  for  multiple  integrals  with  respect  to  a  symmetric  stable  process  (or  measure).  See 
[2]  and  the  references  in  this  paper.  We  consider  multiple  integrals  with  respect  to  a 
completely  asymmetric  stable  process.  We  compute  the  characteristic  function  (cf.).  The  tail 
behaviour  of  the  multiple  integral  will  be  derived  from  the  c.f.  The  behaviour  is  different  from 
the  symmetric  case. 
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LARGE  DEVIATION  RESULTS  FOR  A  CLASS  OF  BANACH  SPACE 
VALUED  MARKOV  CHAINS  WITH  APPLICATION  TO  POPULATION  GENETICS 

G.  J.  Morrow 

University  of  Colorado  at  Colorado  Springs 


The  infinitely  many  alleles  model  in  population  genetics  gives  rise  to  a  probability 
measure  valued  Markov  chain  as  follows.  View  an  allele  as  a  long  sequence  of  nucleotide 
sites.  The  type  t  of  an  allele  will  be  the  site  of  its  most  recent  mutation:  0  £  t  £  1.  Given  a 
population  of  N  alleles  in  generation  n  with  empirical  distribution  v  =  vn  -  £5(t-)/N,  we 
get  a  measure  on  [0,1]  in  generation  n  +  1  defined  by  vn+1  =  £8(tj)/N  where  (tj)  are  i.i.d. 
with  common  distribution  v't  with  an  appropriate  transition  mechanism  v  ->  v'.  We  study  the 


Markov  Chain  (v  )  given  by  v'(A)  =  (1-u) 


(l+s(t))dt/ 


(l+s(t))dt+u  |  A  | ,  A  c  [0,1],  for  a 


mutation  parameter  u  and  a  fitness  parameter  s(t).  Let  these  parameters  admit  the  scaling 
u  =  |ip  and  s(t)  =  J3c(t)  where  P  gives  the  relative  size  of  the  rates  of  mutation  and 
selection.  It  is  assumed  these  evolutionary  effects  are  strong:  N{3  -* «.  The  transition  v-*v' 
has  a  stable  fixed  point  with  density  gg(t)  =  u/[l-Yg(l+s(t))]  when  we  view  v  =  J  <pdv, 
tp  €  A,  as  an  element  of  the  Banach  space  B  =  C(A),  where  A  is  the  compact  subset  of 
C[0,1]  consisting  of  Lipschitz  functions  with  uniform  and  Lipsctiuz  norms  both  bounded  by  1. 
The  logarithm  of  the  expected  time  for  the  chain  [vn)  to  leave  an  open  ball  D  about  gg  in 
B  is  proved  to  be  asymptotically  NPV  for 


V  =  inf  -  2[ 


s(t)(p(t)  -  g0(t))dt  +  | x 


log(p(t)/gg(t))dt]. 


where  the  infimum  is  extended  over  all  densities  p  in  the  complement  of  D.  The  proof  relies 
on  an  infinite  dimensional  extension  of  the  method  of  G.  J.  Morrow  and  S.  Sawyer  appearing 
in  "Large  deviation  results  for  a  class  of  Markov  chains  arising  in  population  genetics",  Ann. 
Prob.  (1989). 
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ASYMPTOTICS  FOR  POINT  PROCESSES 


Gert  Nieuwenhuis 

Katholieke  Universiteit  Brabant 
The  Netherlands 

Let  P  be  the  distribution  of  a  stationary  point  process  on  the  real  line  and  let  P®  be 
its  Palm  distribution.  Our  emphasis  is  on  limit  theorems,  strong  and  weak  (laws  of  large 
numbers,  loglog  laws,  marginal  and  functional  central  limit  theorems).  Two  types  are 
considered,  those  in  terms  of  the  number  of  points  of  the  point  process  in  (0,t]  and  those  in 
terms  of  the  location  of  the  nth  point  right  of  the  origin.  The  former  are  most  easily 
considered  under  P  and  the  latter  under  P®.  General  conditions  arc  presented  that  guarantee 
equivalance  of  either  type  of  limit  theorem  under  both  probability  measures.  As  a  bridge 
between  P^  and  P  we  consider  a  third  probability  measure,  Pp  arising  from  P  by  shifting 
the  origin  to  the  first  point  of  the  process  on  the  right. 

The  obtained  results  for  either  type  of  limit  theorem  are  extended  to  equivalances 
between  the  two  types.  Some  remarks  are  made  on  generalisation  to  non-stationaiy  point 
processes  and  marked  point  processes. 
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LOUD  SHOT  NOISE 


George  L.  O’Brien 
York  University 

(joint  work  with  R.  Doney,  University  of  Manchester) 

Let  (x  }  m<n<ro  denote  the  renewal  epochs  of  a  stationary  renewal  process.  We 
suppose  a  "shot"  occurs  at  each  xn  and  contributes  noise  h(t— xn)  at  time  t,  where 
h:[0,«)  -*  [0,»)  is  a  fixed  function.  We  assume  h(0)  =  1  and  h  is  nonincreasing.  Then  the 
total  noise  at  time  t  is 

5(0  2  “*) 

x 

where  x  runs  through  the  renewal  epochs  up  to  t.  Under  modest  conditions  on  h, 

P[^(t)  >  x  +  1 1  %(t)  >  x]  -+  0  as  x  -*  »,  for  all  t.  Under  more  restrictive  conditions  on  h,  the 
distribution  of  max{l;Crp,...4(T  )}  behaves  in  the  same  way  as  max{ri1,..M7in},  where  the 
Tj/s  are  i.i.d  with  the  same  distribution  as  £(tj). 
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A  PROBLEM  OF  EXPECTATION  AND  VARIANCE 
FOR  A  CLASS  OF  FORMAL  LANGUAGES 

Gabriel  V.  Orman 

Dept,  of  Math.,  University  of  Bra$ov,  2200  Bra§ov,  Romania 


One  of  the  most  interesting  fields  in  which  the  Probability  Theory  can  be  applied  is 
offered  by  the  Mathematical  Theory  of  Formal  Languages.  In  the  last  years  we  have  been 
concerned  with  the  question  of  some  new  characterizations  of  the  process  of  generation  of  the 
words  by  derivations  according  to  formal  grammars.  In  our  studies  the  grammars  constituting 
the  so-called  "Choir  ;<y  hierarchy"  have  been  considered.  To  this  end  we  have  defined  some 
numerical  function:  le  to  characterize  the  derivations  according  to  such  a  grammar  up  to  an 
equivalence  (that  we  call  "derivational  functions"). 

In  our  paper  "Two  problems  concerning  the  probability  of  generation  of  the  words"  (in 
printed),  we  develop  a  particular  study  of  Markov  dependence  in  the  process  of  generation  of 


the  words  by  derivations  according  to  a  phrase-structure  grammar.  Thus,  if  Dx  is  the 
equivalence  class  of  a  derivation  D(x),  x  >  2,  we  denote  by  vx  the  number  of  derivations  of 
D  by  which  a  word  w  is  generated.  Evidently,  w  can  or  cannot  be  generated  into  D  .  If 
it  is,  then  the  probability  that  w  should  be  generated  into  Dx+^  is  denoted  by  y;  if  w  is 


not  generated  into  Dx,  the  probability  that  it  should  be  generated  into  Dx+j 


is  denoted  by 


p.  In  this  way  we  are  in  the  case  when  the  equivalence  classes  D  ,  x  2  2,  are  connected  into 

A 


a  simple  Markov  chain.  In  this  paper  we  continue  this  study  involving  one  of  the  derivational 


functions  namely  the  associated  function  to  a  derivation.  By  means  of  this  function  we 


determine  both  the  probability  that  w  should  be  generated  into  the  class  Dx  and  the 
probability  that  w  should  be  generated  into  the  class  Dj  if  it  was  generated  into  the  class 
D-,i<j.  But  v  is  a  random  variable  that  takes  values  1  and  0  with  probabilities  p  and 
qx  =  l-px  respectively.  We  then  determine  the  expectation  and  the  variance  of  the  random 
variable  giving  the  number  of  derivations  in  n  equivalence  classes  by  which  a  word  w  is 


generated.  Finally,  we  characterize  the  obtained  language  to  the  level  of  this  function. 
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THE  SIZE  OF  DOMAINS  OF  SET-INDEXED  l£VY  PROCESSES. 

Mina  Ossiander 
Oregon  State  University 

This  paper  examines  the  size  of  domains  of  set— indexed  L^vy  processes  with 
independent  increments.  Members  of  this  family  of  infinitely  divisible  random  measures  are 
independent  on  disjoint  sets  and  have  regular  sample  paths  in  the  sense  of  outer  continuity  and 
inner  limits.  An  integral  condition  involving  the  notion  of  majorizing  measure  is  developed 
which  limits  the  size  of  possible  domains.  This  condition  depends  on  both  the  L£vy  measure 
of  the  process  and  the  local  complexity  of  the  domain.  This  work  extends  results  of  Adler  and 
Feigin  (1984)  and  Bass  and  Pyke  (1984).  A  counterexample  of  Adler  and  Feigin  (1984)  is 
generalized  to  show  that  for  p-stable  Levy  processes  the  integral  condition  is  close  to  optimal. 
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ON  THE  TIGHTNESS  OF  THE  PARTIAL  SUMS  OF  A  O-MDONG  SEQUENCE 


Magda  Peligrad 
University  of  Cincinnati 

We  consider  the  partial  sums  of  a  stationary  O-mixing  sequence  of  centered 
L2— integrable  random  variables  having  the  property  that  <  1,  and  we  prove  that  a 
condition  responsible  for  tightness  always  holds.  As  a  consequence,  if  Sn/var(Sn)  converges 
to  a  nondegenerate  distribution,  the  distribution  is  standard  normal. 
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PRECISION  BOUNDS  ON  THE  Lp  NORM  OF  A  QUADRATIC  FORM 
OF  INDEPENDENT  MEAN  ZERO  RANDOM  VARIABLES 


Victor  H.  de  la  Pena 
Columbia  University 

Michael  J.  Klass 

University  of  California,  Berkeley 


Let  Xj,  Xp, ... ,  Xn  be  independent  mean  zero  random  variables, 

S  =  X,  +  X~  +  •  •  •  +  X  .  Fix  p  >  1.  Let  1  £  J  <  n  be  any  index  for  which 
n  1  2  n  r  n  J 

E]Xt  |p  =  max  E|X.|P. 

Jn  1  <,  i<n  1 

Then,  there  exist  constants  0  <  Ap  ,  Bp  <  °» ,  such  that, 

A  E|SnipE[S„-Xj  |PSE|  l  XiXfsB  EISjPEIS^Xj  !p 

n  l^i<j^n  n 

where  Ap  ,  Bp  do  not  depend  on  n  or  the  distributions  of  the  random  variables. 
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ON  ARRIVALS  AND  DEPARTURES 
THAT  SEE  THE  SAME  AVERAGE  IN  BULK  QUEUES 


David  Perry  and  Michael  J.  Magazine 

Department  of  Management  Sciences 
University  of  Waterloo 
Waterloo,  Ontario 
Canada  N2L  3G1 


We  consider  a  general  multi-server  queueing  system  with  batch  arrivals  in  which  the 
arrival  process  is  a  renewal  process.  For  this  general  system  we  show  that  there  are  special 
relationships  between  the  limiting  probability  law  at  departure  times  and  the  arrival  times.  For 
the  special  cases  of  Poisson  arrivals  or  exponential  service  time,  we  generalize  these 
relationships  to  the  stationary  distribution  in  time  average. 
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NON-LINEAR  FILTERING  WITH  SMALL  NOISE  COEFFICIENTS 

Jean  Picard 

INRIA-Sophia  Antipolis,  2004  route  des  Lucioies,  06560  Vaibonne,  FRANCE 

Consider  the  non-linear  filtering  problem  where  the  signal  Xt  and  the  observation  Y 
are  solutions  of 

'  dXt  =  P(t,Xt)dt  +  o  ( t,Xt)dWt  +  v€  7(t>Xt)dBtl 
.  dYt  =  h(t,Xt)dt  +  y€  dBt 

with  two  independent  Brownian  motions  W  and  Bt.  The  aim  of  this  work  is  to  study 
asymptotically  the  conditional  law  of  X  given  Y  as  e  -*  0.  Firstly,  we  find  conditions  under 
which  the  variance  of  the  conditional  law  is  of  order  e  for  sufficiently  large  times;  these 
conditions  include  the  cast  where  x  w  h(t,x)  is  one— to— one  (the  easiest  case)  but  also  some 
other  situations;  this  problem  is  linked  with  the  existence  of  non-linear  observers  for  the 
associated  deterministic  system  (obtained  for  e  =  0)  and  use  the  notion  of  detectability. 
Secondly,  we  check  that  under  some  additional  assumptions,  the  conditional  law  is 
asymptotically  Gaussian;  in  particular  this  provides  a  rigorous  justification  of  the  famous 
extended  Kalman  filter,  this  suboptimal  filter,  based  on  a  linearization  of  the  system,  is  proved 
nearly  optimal  in  this  case;  however,  we  also  describe  other  cases  where  it  behaves  badly. 

By  a  change  of  time  scale,  we  can  also  study  the  situation  where  the  orders  of 
magnitude  of  the  signal  and  observation  noises  are  different;  with  this  remark,  we  can  recover 
as  a  particular  case  th<*  results  of  previous  works  such  as  [1,2].  The  mathematical  tools  which 
are  needed  include  the  stability  of  stochastic  systems,  some  estimates  of  large  deviations,  some 
changes  of  probability  measures  and  the  stochastic  calculus  of  variations. 
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CENTRAL  LIMIT  AND  RECURRENT  BEHAVIOR  OF  DIFFUSIONS 
IN  A  RANDOM  ENVIRONMENT 


Ross  Pinsky 
Technion 

and 

Mark  Pinsky 
Northwestern  University 


Let  (X(t),c(t))  be  the  Markov  process  on  the  state  space  IR^  x  {1,...,N}  with  the 
infinitesimal  generator  L  =  +  G  +  V(x,c)  where  G  is  the  infinitesimal  matrix  of  a 

Markov  chain  on  {1,...,N},  A  is  the  Laplacian  of  IR^  and  V(x,a)  is  a  homogeneous  vector 
field  with  mean  zero  (w.r.t.  the  invariant  measure  of  (o(t):  £0}). 

Theorem  1:  When  e  [  0,  e  X(t/ez)  converges  weakly  to  a  diffusion  process  on  Ra  with 
generator  Lg  =  ^A  + 

I 

Markov  chain. 


I 


c~oV(x,a)V(x,(3)  where 


'ap 


cap  is  the  recurrent  potential  kernel  of  the 


Theorem  2:  If  d  =  2  the  process  (X(t):  t>0}  is  recurrent. 

Theorem  3:  If  d  =  3  the  process  { X(t):  £0}  is  transient  if  sup|V(x,c)|  is  sufficiendy 
small. 
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A  PROBABILISTIC  ANALYSIS  OF  UPSCHTTZ  POINTS 
FOR  NON  DIFFERENTIABLE  A*  FUNCTIONS 

J.  M  Anderson  and  L.  D.  Pitt 
University  of  Virginia 


* 

A  continuous  periodic  function  f(x)  is  in  A  if 


sup  (1/h)  |f(x+h)  +  f(x-h)  -  2f(x)  |  <  «. 
h>0,x 

* 

A  functions  may  be  nowhere  differentiable,  as  in  the  case  with  the  Weierstrass  function 

f(x)  :=  ^2_ncos  2nx, 

* 

but  must  satisfy  Lipschitz  conditions  on  a  large  set.  For  f  e  A  we  use  the  martingale  of 
dyadic  difference  quotients 


Dnf(x)  :=2n^[f((k+l)/2Vf(k/2n)]  •  1 


r  (x), 
n,k 


where  I  >  =  [k/2n,  (k+lj/2n),  to  analyze  the  sets 

n,K 


A(c)  =  {x:  lim  sup|(l/n)(f(x+h)-f(x))|  £c} 
h-0 


B(c)  =  {x:  lim  sup  | 
h-0 


Upper  and  lower  bounds  for  the  size  of  the  sets  A(c)  and  B(c)  are  obtained,  which  in  some 
cases  give  exact  measure  functions. 
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ESTIMATION  OF  RELIABILITY  OF  STRESS-STRENGTH  MODELS 
BASED  ON  SAMPLE  PATHS 


Nader  Ebrahimi  and  T.  Ramalingam 

Northern  Illinois  University 
Division  of  Statistics 
Department  of  Mathematical  Sciences 
DeKalb,  IL  60115 

The  stress  and  strength  patterns  of  many  stochastic  systems  are  time  dependent. 
Consequently,  traditional  non— dynamic  models  for  stress  and  strength  are  inadequate  for 
assessing  the  reliability  of  such  systems.  Furthermore,  it  is  more  effective  to  monitor  systems 
over  time  and  estimate  the  reliability  or  predict  the  time  to  failure  of  the  system  in  advance  of 
failure,  than  assuming  that  failure  times  are  available  for  inference  about  the  system. 

In  this  paper,  we  propose  a  general  stochastic  model  for  stress  and  strength  of  systems 
and  illustrate  related  inferences  in  some  special  cases  including  the  Brownian  Motion  Model. 
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LIGHT  AND  HEAVY  TRAFFIC  LIMITS  FOR  TWO 
FINITE  CAPACITY  QUEUEING  SYSTEMS 


Martin  L  Reiman 

AT&T  Bell  Laboratories 
Murray  Hill,  NJ  07974 

Burton  Simon 

University  of  Colorado  at  Denver 
Denver,  CO  80202 

We  consider  two  related  queueing  systems  which  have  two  classes  of  customers: 
infinite  capacity  (IC)  and  finite  capacity  (FC).  Both  customer  classes  have  Poisson  arrival 
processes.  Upon  their  arrival,  IC  customers  enter  the  main  queue,  which  is  served  by  a  single 
server.  There  is  a  pool  of  N  tokens  associated  with  FC  customers,  and  an  FC  customer 
needs  a  token  to  enter  the  main  queue.  Customers  in  the  main  queue  are  served  in 
first-come-first-served  order.  Service  time  distributions  are  general  and  (possibly)  different 
for  each  class.  In  Model  1,  if  no  token  is  available  an  arriving  FC  customer  waits  in  a 
secondary  queue  until  a  token  becomes  available.  In  Model  2,  if  an  arriving  FC  customer 
finds  no  token  available  it  leaves  the  system. 

We  determine  the  asymptotic  behavior  of  this  system  in  light  and  heavy  traffic.  Our 
light  traffic  asymptotics  consist  of  the  first  N  +  1  derivatives  of  sojourn  time  moments  with 
respect  to  the  arrival  rate,  at  an  arrival  rate  of  zero.  The  heavy  traffic  limit  of  IC  customers  is 
the  ame  in  both  models:  a  one  dimensional  diffusion  similar  to  a  Bessel  process.  The  heavy 
traffic  limit  of  IC  customers  is  the  same  in  both  models:  a  one  dimensional  diffusion  similar  to 
a  Bessel  process.  The  heavy  traffic  limit  of  FC  customers  in  Model  1  '  a  one  dimensional 
reflected  Brownian  motion.  These  asymptotics  can  be  combined,  by  interpolating  between 
them,  to  yield  an  approximation  for  all  arrival  rates. 
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APPROXIMATIONS  FOR  DIFFUSION  PROCESS  BOUNDARY  HITTING  TIMES 


Gareth  Roberts 

University  of  Nottingham,  England 

Methods  for  approximating  distributions  of  boundary  hitting  times  for  diffusion 
processes  will  be  considered.  First,  limit  theorems  for  the  distribution  of  a  diffusion 
conditioned  not  to  hit  a  time-invariant  boundary  are  established,  leading  to  tail  estimates  for 
the  hitting  time  of  a  constant  boundary.  Secondly,  these  results  are  used  to  give 
asymptotically  correct  estimates  of  time  dependent  boundary  hitting  times. 

Applications  of  these  results  in  statistics  will  be  briefly  discussed. 
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BROWNIAN  MOTION  IN  A  WEDGE  WITH  VARIABLE  SKEW  REFLECTION 

L.  C  G.  Rogers 
University  of  Cambridge 

In  recent  years,  there  has  been  considerable  interest  in  the  behaviour  of  Brownian 
motion  in  a  wedge  with  skew  reflection  on  the  two  sides  of  the  wedge,  the  direction  of 
reflection  being  assumed  constant  on  the  two  sides  of  the  wedge.  This  talk  will  discuss  the 
reachability  of  the  vertex  when  the  direction  of  reflection  is  allowed  to  vary  along  the  edges. 
By  techniques  of  complex  analysis,  it  is  possible  to  find  a  necessary  condition,  and  a  sufficient 
condition,  for  the  vertex  to  be  reachable.  These  conditions  are  of  a  geometric  nature;  the  gap 
between  them  is  evidently  quite  small! 
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A  REPRESENTATION-INDEPENDENT  INTEGRAL 
FOR  SET-VALUED  RANDOM  VARIABLES 

David  Ross 

Dept,  of  Math,  and  Stat. 

University  of  Minnesota 
Duluth,  MN  55812 

If  %  is  an  infinite-dimensional  Banach  space,  and  F  and  F  are  random  subsets  of 
X  (i.e.,  random  variables  taking  values  in  7Xx))»  then  it  can  happen  that  F  and  F  have  the 
same  distribution  hut  different  integrals  (in  the  sense  of  the  Aumann  integral).  Debreu's  notion 
of  integration  addresses  this  problem,  at  the  cost  of  restricting  the  class  f  integrable  sets. 

In  this  talk  I  discuss  several  natural  definitions  for  the  integral  of  a  general  set— valued 
random  variable  which  depend  solely  on  the  r.v.'s  distribution.  Remarkably,  these  definitions 
all  produce  the  same  integral. 
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STOCHASTIC  GEOMETRY  AND  RELATIVISTIC  PATH  INTEGRALS 

Sisir  Roy 

Indian  Statistical  Institute 
Calcutta  700035,  INDIA 

Feynman  examined  the  Klein— Gordon  and  Dirac  equations  by  the  method  of 
introducing  a  fifth  coordinate  u  which  played  the  role  of  a  proper  time.  The  jet  extension  of 
the  Finslerian  fields  displays  the  idea  of  covariance  of  the  physical  field  theory  with  respect  to 
general  changes  of  the  parameter  defined  along  the  trajectories.  It  seems  that  this  offers  an 
attractive  and  quite  obvious  opportunity  of  consistent  construction  of  dynamics  of  scalar, 
spinor,  gauge  etc.  physical  fields  with  respect  to  the  parameter  l  We  propose  to  call  and 
consider  such  parameter  t  as  the  evolution  time  parameter.  A  field  can  be  said  to  be 
evolutionary-dynamic  if  it  is  given  by  a  Lagrangian  involving  the  first  derivative  of  the  field 
with  respect  to  the  evolutionary  time  parameter.  This  velocity  can  take  arbitrary  values  at  a 
certain  point  with  fixed  x1  and  t  in  the  four  dimensional  differentiaDle  Finslerian  Manifold. 
Then  the  metric  tensor  will  depend  on  the  random  velocity  variable.  Now  by  using  such 
covariant  parametrization  and  arbitrary  velocity  vector,  the  description  of  Feynman's  path  as 
trajectories  is  possible  even  for  an  individual  particle. 
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ON  THE  FINITENESS  OF  THE  MOMENTS  OF  FIRST-CROSSING  TIMES 
FROM  NONLINEAR  RENEWAL  THEORY 


Valeri  T.  Stefanov 


Bulgarian  Academy  of  Sciences 
Sofia,  Bulgaria 
and 

Odense  University 
Odense,  Denmark 


The  following  stopping  time  covers  a  number  of  important  special  cases  occurring  in 
nonlinear  renewal  theory  and  sequential  analysis: 

x  =  inf{n:  ng(n_1Zn)  t  a},  a  >  0  , 

where  Zfl  is  a  random  walk  with  a  positive  drift  (p.)  and  g(  • )  is  a  positive  function  which 
is  twice  continuously  differentiable  in  a  neighbourhood  of  p.  We  show  that  all  moments  of  t 
are  finite  in  the  case  when  the  random  walk  follows  a  distribution  from  any  noncurved 
exponential  family.  This  is  also  true  for  the  continuous-time  analogue  of  T. 


88 


PRODUCT  FORM  IN  MIGRATION  PROCESSES 
WITH  BATCH  MOVEMENTS 

W.  Henderson 
University  of  Adelaide 

P.  Taylor 

University  of  Western  Australia 

Walrand  [6]  and  Pujolle  [4],  [5]  have  recently  derived  the  equilibrium  colony  size 
distribution  for  certain  migration  processes  in  which  customers  can  arrive  and  depart  from 
colonies  in  batches.  In  this  talk  we  will  define  a  general  class  of  migration  processes  which 
includes  the  models  of  both  Walrand  and  Pujolle  and  discuss  a  method  by  which  the 
equilibrium  colony  size  distribution  may  be  found. 

It  turns  out  that  for  processes  satisfying  appropriate  conditions  the  equilibrium 
distribution  has  a  form  which  is  analogous  to  the  product  form  which  has  been  observed  in 
migration  processes  with  single  customer  movements  by,  among  others,  Jackson  [2],  Baskett, 
Chandy,  Muntz  and  Palacios  [1]  and  Kelly  [3]. 
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SUMS  AND  MAXIMA  OF  STATIONARY  SEQUENCES 


K.  F.  Turkman 
CEAUL 

University  of  Lisbon 


Let  (Xn)  be  a  stationary  uniform  mixing  sequence  with  distribution  function  F(.) 


and  mixing  coefficient  <p(j)  satisfying 


Let  {Xn)  be  the  associated  i.i.d.  sequence  with  the  same  distribution  function.  Let  (an.bn) 

-1  r11  *  -1  * 

and  (cn,dn)  be  constants  such  that  ^  (Xj-b^  ^  cn  (maxi<i<nx;~^n)  converge 

in  distrubution  to  S  and  M  respectively,  where  S  is  stable  with  index  0  <  a  £  2  and  M 
is  a  non  -degenerate  random  variable.  In  this  paper,  under  a  local  dependence  condition  and 
(1),  thit  joint  limiting  distribution  of 

I  cn1(,S,W 

is  studied  and  the  possible  limiting  forms  are  given.  It  is  shown  that,  except  for  the  case  when 
S  is  stable  with  index  a  <  2  and  M  is  max-stable  with  the  same  index,  the  normalized 
sums  and  maxima  are  independent. 
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NONPARAMETRIC  INFERENCE  FOR  A  DOUBLY  STOCHASTIC  POISSON  PROCESS 


Klaus  J.  Udkal 

Department  of  Statistics 
University  of  Kentucky 
Lexington,  Ky  40506 


Consider  a  doubly  stochastic  Poisson  process  whose  intensity  is  given  by 

Xt  =  a(Zt),  where  a  is  an  unknown  nonrandom  function  of  another  (covariate)  process  Z 

Only  one  continuous  time-observation  of  counting  and  covariate  process  is  available.  The 
z 


function  A(z)  = 


a(x)dx  is  estimated  and  the  normalized  estimator  is  shown  to  converge 


6 


weakly  to  a  Gaussian  process  as  time  approaches  infinity.  Confidence  bands  for  A  are  given. 
Tests  for  independence  from  the  covariate  process  are  proposed. 


Key  words  and  phrases:  Optical  Communications,  Doubly  Stochastic  Poisson  Processes, 
Nonparametric  Functional  Estimation,  Method  of  Sieves,  Martingale  Central  Limit  Theorem, 
Semimartingales,  Gaussian  Processes. 
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SOME  PROBABILISTIC  PROBLEMS  ARISING  FROM 
AN  OPTIMAL-PATH  ALGORITHM 


J.  van  den  Berg 

Centre  for  Mathematics  and  Computer  Science 
Kmislaan  413 
P.O.  Box  4079 
1009  AB  Amsterdam 

Consider  a  finite  graph  G.  Suppose  that  with  each  edge  e  two  types  of  cost,  c(e) 
and  d(e)  are  associated. 

Let  A  and  B  be  vertices  of  G.  The  path  from  A  to  B  which  minimizes  cost  c 
can  be  found  by  Dykstra's  algorithm.  The  running  time  of  this  algorithm  is  polynomially 
bounded  in  the  number  of  edges  of  G.  If  one  wants  to  minimize  cost  c  under  an  additional 
constraint  on  cost  d,  or,  more  generally,  if  one  wants  to  minimize  an  increasing  but 
non-linear  function  of  the  total  cost  c  and  d,  this  can  be  done  by  a  generalised 
(multiple-labelling)  Dykstra  algorithm.  The  running  time  is  now  closely  related  to  the  number 
of  efficient  paths  between  two  given  vertices  (a  path  is  efficient  if  there  is  no  other  path  of 
which  both  cost  c  and  cost  d  are  smaller).  It  is  not  difficult  to  construct  an  example  where 
this  grows  exponentially  with  the  number  of  edges.  We  take  the  costs  d(e)  and  c(e) 
stochastic  and  study  the  average  number  of  efficient  paths.  This  leads  to  more  general 
interesting  probabilistic  problems  of  a  combinatorial  nature. 
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RANDOM  WALK  APPROXIMATIONS  OF  SKEW  AND  STICKY  BROWNIAN  MOTION 

J.  A.  M.  Van  Per  Weide 
Delft  University  of  Technology 


Let  (Sm)m>f)  be  a  Markov  chain  on  Z  which  behaves  outside  0  like  symmetric 
random  walk.  From  0  the  chain  makes  a  transition  to  a  state  k  according  to  a  probability 
distribution  (p^)^  j  on  2  with  finite  first  moment.  Consider  the  process  Xn  =  (Xn(t))^Q, 
n  =  1,2,...  defined  by: 


Xn<‘>  = 

Harrison  and  Shepp  [1]  claim  weak  convergence  of  the  sequence  Xn  to  skew  Brownian 
motion  with  a  parameter  depending  on  the  distribution  '^c  ’*^1  present  a  proof  of 

this  result  using  semigroup  methods.  As  a  special  case  of  the  more  general  problem  where  the 
probability  distribution  (Pjc)jc6  j  depends  on  n,  we  will  discuss  a  random  walk  approximation 
of  sticky  Brownian  motion. 
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NEW  RESULTS  FOR  OPTIMAL  SWITCHING 
BETWEEN  A  PAIR  OF  BROWNIAN  MOTIONS 

Robert  J.  Vanderbei 

AT&T  Bell  Laboratories 
Murray  Hill,  NJ  07974 

* 

We  consider  the  problem  to  find  the  T  that  achieves  the  following  sup 

v(x)  =  suPt  Ex  f(BT1(x)  *BT2(t))  * 

1  2 

where  B(s)  =  (B„  ,B„  )  is  a  two-parameter  Brownian  motion,  f  is  a  pay-off  function 
S1  s2 

2 

defined  on  the  boundary  of  the  unit  square  in  R  ,  x  is  the  first  exit  time  from  the  unit  square 
and  T  =  (T(t)  =  (TjCO.TjCt)):  0  <  t  <  «}  is  an  optional  increasing  path. 

In  this  talk  we  quickly  review  old  results  and  then  discuss  new  ones. 

This  is  joint  work  with  Avi  Mandelbaum  and  Larry  Shepp. 
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LOGLOG  LAWS  FOR  GEOMETRIC  AND  FULL  SEQUENCES 


Wim  Vervaat 

Catholic  University 
Nijmegen,  The  Netherlands 

Why  are  functional  loglog  laws  like  Strassen's  proved  first  for  geometric  subsequences? 
Is  the  subsequence  result  indeed  weaker?  The  language  of  dynamical  systems  gives  the  best 
insight. 

Let  a  group  G  ict  continuously  on  a  metric  space  S.  Let  y  be  a  measurable  function 
from  IR  into  G  whicn  is  regularly  varying,  i.e.,  limt_^y(t+a)(7(t))-*  *  5(a)  in  G 
(consequently,  5(a+b)  =  8(a)5(b),  a  h  5(a)  is  continuous  and  the  convergence  is  locally 
uniform).  Let  f  e  S  be  fixed  and  such  that  for  each  sequence  t  -» »  in  !R  the  sequence 
7(t  )f  in  S  has  a  convergent  subsequence,  and  let  K  be  the  set  of  limit  points  in  S  of  y(t)f 
as  t  -*  «.  For  real  a  >  0,  let  K  cK  be  the  set  of  limit  points  of  y(na)f  as  n  neW. 

Questions:  when  do  we  have  K  =  K;  for  how  many  a  can  we  have  K  *  K?  Here 

a  a 

is  the  translation  of  Strassen's  loglog  law  for  Brownian  motion  B  into  this  language: 

S  =  C[0,«) , 

f  =  B(  •  ,to)  for  a  fixed  sample  point  to  , 

tflog  t)f  =  -~==z> 
y/2l  loglogt 

5(log  a)f  =  a_1/2f(a-) , 

00 

K  =  {f  6  S:  f(0)  =  0,  f  absolutely  continuous,  (f(t))2dt  SI). 

0 

Several  other  functional  loglog  laws  to  which  this  can  be  applied  are  being  explored  now. 
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EXPECTED  AMOUNT  OF  OVERFLOW  IN  A  FINITE  DAM 


P.  R.  Vital 

Department  of  Mathematics 
Western  Carolina  University 
Cullowhee,  NC  28723 

In  a  finite  dam  model  with  random  inputs  and  random  outputs  the  expected  amount  of 
overflow  in  a  given  time  is  analyzed  using  imbedding  method.  A  second  order  differential 
equation  for  the  expected  amount  of  overflow  is  obtained  and  its  complete  solution  is  also 
derived.  When  the  input  is  random  and  the  output  is  proportional  to  the  content  of  the  dam,  a 
second  order  differential  equa  on  of  Kummer's  type  is  obtained  and  its  closed  form  solution  is 
derived  in  terms  of  confluent  hypergeometric  functions.  In  addition  to  random  outputs,  if  a 
deterministic  exponential  release  policy  is  also  followed  we  arrive  at  a  third  order  differential 
equation  for  the  expected  amount  of  overflow  and  its  solution  is  in  terms  of  beta  functions  and 
degenerate  hypergeometric  functions  of  two  variables. 
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THE  SUBSTITUTION  METHOD  FOR  BOND  PERCOLATION 
CRITICAL  PROBABILITY  BOUNDS 


John  C  Wiennan 

Mathematical  Sciences  Department 
Johns  Hopkins  University 
Baltimore,  Maryland 


A  "substitution  method"  for  deriving  bounds  for  critical  probabilities  of  bond 
percolation  models  will  be  presented.  In  this  method,  a  lattice  is  compared  with  another  lattice 
for  which  the  critical  probability  is  exactly  known.  Portions  of  one  lattice  are  substituted  into 
the  other,  while  an  appropriate  transformation  is  applied  to  the  parameters  in  the  bond 
percolation  model.  The  transformation  is  determined  from  a  stochastic  ordering  on  probability 
measures  derived  from  finite  random  graphs  corresponding  to  portions  of  the  lattices.  The 
method  theoretically  provides  sequences  of  upper  and  lower  bounds  for  the  bond  percolation 
critical  probability.  Both  upper  and  lower  bounds  for  the  bond  percolation  critical  probability 
of  the  Kagome  lattice  are  substantially  improved  by  an  application  of  this  method. 
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ON  THE  QUASIREVERSIBILITY  OF  A  MULTICLASS 
BROWNIAN  SERVICE  STATION 

J.  M.  Harrison 
Stanford  University 

R.  J.  Williams 

University  of  California  at  San  Diego 

The  object  of  study  in  this  paper  is  a  Brownian  model  of  a  multiclass  service  station. 
Such  Brownian  models  arise  as  heavy  traffic  limits  of  conventional  queueing  models  in  which 
several  different  types  or  classes  of  customers  are  processed  through  a  common  service 
facility.  Assuming  that  the  Brownian  service  station  is  initialized  with  its  stationary 
distribution,  four  different  model  characteristics  are  shown  to  be  equivalent,  and  the  station  is 
said  to  be  quasireversible  if  those  equivalent  conditions  pertain.  Three  of  the  four  conditions 
characterize  the  vector  output  process  from  the  Brownian  service  station,  and  our  definition  of 
quasireversibility  parallels  precisely  that  proposed  by  F.  P.  Kelly  for  conventional  queueing 
models.  The  last  of  our  four  conditions  is  expressed  directly  in  terms  of  primitive  station 
parameters,  so  one  may  easily  determine  from  basic  data  whether  or  not  a  Brownian  station 
model  is  quasireversible.  In  future  work  we  hope  to  show  that  a  network  of  quasireversible 
Brownian  service  stations  has  a  product  form  stationary  distribution. 
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EQUIVALENT  MARTINGALE  MEASURES  AND  NO-ARBITRAGE 
IN  STOCHASTIC  SECURITIES  MARKET  MODELS 


Robert  G.  Dalang 
Department  of  Statistics 
University  of  California  at  Berkeley 
Berkeley,  CA  94720 

Andrew  Morton 

Dept,  of  Information  and  Decision  Sciences 
University  of  Illinois  at  Chicago 
Chicago,  IL  60680 

Walter  Willinger 
Bellcore 

Morristown,  NJ  07960 


We  characterize  those  vector-valued  stochastic  processes  (with  a  finite  index  set  and 
defined  on  an  arbitrary  stochastic  base)  which  can  become  a  martingale  under  an  equivalent 
change  of  measure. 

This  question  is  important  in  a  widely  studied  problem  which  arises  in  the  theory  of 
finite  period  securities  markets  with  one  riskless  bond  and  a  finite  number  of  risky  stocks.  In 
this  setting,  our  characterization  gives  a  criterion  for  recognizing  when  a  securities  market 
model  allows  for  no  arbitrage  opportunities  ("free  lunches").  Intuitively,  this  can  be 
interpreted  as  say  "if  one  cannot  win  betting  on  a  process,  then  it  must  be  a  martingale  under 
an  equivalent  measure",  and  provides  a  converse  to  the  classical  notion  that  "one  cannot  win 
betting  on  a  martingale". 
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SOME  NEW  VAPNDC  -  CHERVONENK3S  CLASSES  OF  SETS 


G.  Stengle  and  J.  E.  Yukich 

Department  of  Mathematics 
Lehigh  University 
Bethlehem,  Pennsylvania  18015 


This  paper  provides  some  new  examples  of  Vapnik  -  Chervonenkis  (VC)  classes  of 
sets  and  shows  that  these  classes  fit  naturally  within  the  framework  of  algebraic  geometry.  In 
particular,  semialgebraic  sets  are  linked  to  VC  classes  and  are  used  to  manufacture  relatively 
large  VC  families  of  of  positivity  sets,  thus  extending  a  result  of  Dudley  (Ann.  Prob.  6 
(1978):  899-929).  The  main  tool  here  is  the  quantifier  elimination  theorem  of 
Tarski-Seidenberg.  Similarly,  tools  from  the  theory  of  analytic  functions  are  used  to  show 
that  certain  special  analytic  families  of  positivity  sets  are  VC. 


